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NOUVELLES MÉTHODES 


POUR LA DÉTERMINATION 


DE L’ORBITE DES COMÈTES. 


SECOND SUPPLÉMENT. ' *' 

X,( ■ J 


D >y ^ ^ 

ans le supplément qui précède , nous avons pev^âCtsVfiné et 
simplifié la méthode exposée dans notre premier Mémoire, de 
manière qu’elle devra donner dans tous les cas une approximation 
suffisante pour la pratique, surtout si on l’applique à des obser- 
vations effectives et qui n’aient point été soumises à l’interpolation. 
Cette méthode serait susceptible d’un plus grand degré de pré- 
cision, si on voulait admettre les quantités du quatrième ordre 
dans le développement des-coordonnécs , ce qui rendrait les 
cqualious (e) et (g), (pag. 5 et 6), exactes jusqu’aux quantités 
du second ordre inclusivement ; mais l’introduction de ces termes 
augmenterait beaucoup la longueur des calculs, et rendrait les 
formules moins élégantes. 

Nous nous proposons maintenant d’exposer une seconde mé- 
thode par laquelle on peut assez facilement obtenir la solution 
la plus exacte dont le problème est susceptible , et qui a l’avan- 
tage de se rapprocher de plusieurs méthodes connues, telles que 
la méthode de M. Olbers et celle que M. Ivory a donnée dans 
les Transactions philosophiques de la Société royale de Londres 
an i8i4, première partie. 

Nous considérons d’abord le secteur décrit par le rayon vec- 
teur delà comète, dans le temps écoulé de la première à la Iroi- 
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sième observation , et nous cherchons dans quel rapport la corde 
de ce secteur est divisée par le rayon vecteur qui répond à la 
seconde observation. Oiï sait, d’apAs Newton, que ce rapport 
diffère en général très peu du rapport des temps écoulés dans les 
intervalles des observations, et les approximations données par 
MM. Olbers et Ivory supposent que la différence entre ces deux 
rapports peut toujours être négligée. 

Mais cette supposition n’est pas suffisamment exacte dans plu- 
sieurs circonstances. Aussi la rectification successive de la valeur 
attribuée au rapport des deux parties de la corde, caractérise la 
différence essentielle entre la méthode nouvelle que nous pro- 
posons , et celles qu’ont données MM. Olbers et Ivory. 

En général, nous distinguons, dans cette seconde méthode, 
trois degrés successifs d'approximation ; le premier , dans lequel 
on néglige lcs_ {ermes„du second ordre en 0 et 0', ne donne 
que des résultats peu exacts, mais qui suffisent pour diriger les 
hypothèses dans les calculs suivans. 

Le second degré d’approximation , dans lequel on a égard aux 
quantités du second ordre, est à peu près le même que celui 
qu’on obtient par notre première méthode ; il revient encore à 
celui que donneiu W méthodes de MM. Olbers et Ivory , dan* 
tous les cas où l’on peut supposer amis erreur sensible que la 
corde est partagée en raison des temps. 

Enfin le troisième degré d’approximation est celui où l’on 
tient compte des quantités du troisième ordre; il s’obtient très 
facilement au moyen de l’approximation précédente , et l’orbite 
qui en résulte satisfait eu général aussi exactement qu'il est pos- 
sible aux observations données. 

Cette méthode cependant ne s’appliquerait plus avec le même 
succès, si l’orhitc réelle s’écartait sensiblement de la figure pa- 
rabolique; c’est ce qu’on reconnaîtrait par les éiémens trouvés 
de l’orbite parabolique, lesquels satisferaient toujours exactement 
aux deux observations extrêmes, mais ne représenteraient l’ob- 
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servation moyenne qu’avec des erreurs plus ou moins considé- 
rables. Dans ce cas , il faudrait avoir recours aux me'lbodes con- 
nues des astronomes pour la détermination des orbites elliptiques 
ou hyperboliques. D’ailleurs, si on n’a besoin que d’une première 
approximation , on pourra en trouver une par notre première 
méthode, en n’omettant aucune des équations du problème, et 
déterminant le demi-grand axe n par la dernière de ces équa- 
tions, savoir ~ ^ = »»'• -f- «’• -f- p‘. 

Pour rendre ces recherches plus complètes , j’ai cru qu’il ne 
serait pas inutile d’exposer sommairement une troisième mé- 
thode, où l’on suppose connu, pour un instant donné, la lon- 
gitude et la latitude apparentes de la comète, avec leurs cocf- 
ficiens différentiels du premier et du second ordre, pris par 
rapport au tems. On a reproché à cette Qiélhode, la longueur 
des calculs d’interpolation nécessaires pour déduire ces coefü- 
ciens d’un certain nombre d’observations données ; nous avons 
fait remarquer de plus, page iv du premier Mémoire, que les 
résultats de l’interpolation peuvent devenir d’autant plus inexacts, 
qu’on combine ensemble un plus grand nombre d’observations. 
Quoi qu’il en soit, toutes les fois que les coefliciens différentiels 
pourront être déterminés avec une précision suffisante, celle 
méthode sera la plu» simple et la plus exacte de toutes. C’est 
pourquoi nous la donnons ici avec des explications qui en faci- 
literont l’usage et la feront mieux apprécier. 

Je n’ajouterai rien aux méthodes que j'ai déjà données dans 
les deux Mémoires précédées et surtout dans le second, pour 
trouver l’orbite corrigée qui satisfait, le plus exactement qu’il 
est possible, à un nombre quelconque d’observations, faites dans 
un espace de temps plus ou moins long. 11 me semble que jusqu'à 
la publication de ccs Mémoires, qui a eu lieu en i8o5el j8o6, 
la question n’avait été traitée que d'une manière imparfaite et 
par de simples tâtonnemeus. J’ai donné le premier deux méthodes 
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sures pour obtenir la solution à la fois la plus simple et la pins 
exacte, savoir : la méthode des corrections indéterminées qui n’avait 
encore reçu qu'un petit nombre d’applications, et la méthode 
des moindres carrés , qui paraissait alors pour la première fois. 


SECONDE MÉTHODE. 

• 

i. Soit S le centre du Soleil (Gg. i ), Sy, l'axe des x dirigé vers 
le premier point d 'Âries; T°, T, T' les trois lieux de la Terre 
correspondans aux temps j — 8 , /,/- f-8'; soient C“, C, C', les 
lieux observés de la comète dans les trois mêmes ins tans, et 
enGn Q% Q, Q', les projections de ces points sur le plan de 
l’écliptique. Nous appellerons : 

!A*, A, A', les trois^ longitudes héliocentriques de la Terre; 

R%R, R-, -des trois rayons vecteurs ST", ST, ST'; 

a*, a, a', les trois longitudes géocenlriques de la comète; 

r *, r, d, les trois rayons vecteurs SC*, SC, SC'; 

b’, b, b' y les trois latitudes géocenlriques de la comète; 

<r°, v, <t\ les trois distances accourcies TQ% TQ, TQ'; 

c, la corde C*C de l’arc parcouru pendant le temps ô-f-S'. 

Le plan de l’orbite passant toujours par le Soleil, le rayon vecteur 
CS coupera la corde C*C' en un point D qui aura pour projection 
le point E. 

On sait que le rapport des deux parties de la corde C*D, DC', 
est à très peu près égal à celui des temps 8, 8' ; nous supposerons 
que ce rapport est exactement celui de A à A', et nous prendrons 
A-f-A' = i , ce qui donnera C*D=Ae et C'D=A'c. 

Enfin nous supposerons CD=er, ou SD=(i — e)r, e devant 
être une quantité très petite du second ordre, par rapport à 8 
et 8' qui sont censés du premier ordre. 

a. Pour trouver le rapport de A à A', qui sera le fondement de 
notre analyse, nous supposerons d’abord l’orbite parabolique. 
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Soit la distance périhélie SP=n ( fig. a ), l'anomalie au point 
C on l'angle PSC=4<, l’angle C°SC=£, l’angle CSC' = £'; 
•oit de plus tang-j-vf ,— u > tan gï^=?'» tangi£' = j/; on aura par 
les propriétés de la parabole , 


r =n(i+ ü *), 

r*=n(x-f u»»), 
/ / = n(i+u'*). 


u 1 >=taDg(i4— JC) 

u'=tang(i+ + iî') : 


a — y 

x *4- yu 1 
u+y 

I — y'u’ 


... rÇ. -1- •/) 
(i +y«)“ 
._ >•(■ +y'*> 
C-yu)*’ 


Dans les triangles C°DS, C'DS, on a 


donc 


C S D = 


r° >in C 
sinD * 


C'D = 


r’sin f 
sinD * 


^ r° si n C y ( i — y'u)* 

V — ~ V liiiS' ÿ 7 ’ (i -J- y u)“" 


Dans l’orbite elliptique, on a le quarré de la vitesse 

^ = 2 — i; dans l’orbite parabolique où a est infini, on a 

simplement , ou rfr * — = y Substituant les valeurs...- 
d-\, = ^-~r;,djr= ail «du, on aura <//-*-f*r*<f4‘ ==: 4n*( I + u')du* 
c=jdt‘. De là résulte = du { i -4- u*) , et en intégrant à 

compter du périhélie, 

= « + ï“ s î 

_ jSÜfPIRI. 7 •• -;v 

mettant / -f- 6' an lieu de t et u' au lieu de u, on aura 

ou en substituant la valeur de u r , 

< _ y'0+u > )V . . , y’ 1 (i-4-ü>A , 

ÿ'(an') (x — y'u)* \ ' î — y u / 

Cette équation est relative au point C' ; on aurait semblablement 
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pour le point C* 

fl __ vC» 4- ’t'Y / , , »*(i + u ')\ 

|/(ariO (i -f-yu)* \ ' i + yj- )' 

Soit pour un moment Z = ■ ' . Z ' == — — ; on aura les 

r i-yu' 

deux équations • 


= + i>^- 

De là on déduit, en rejetant seulement les quantités du cin- 
quième ordre, 

= Z'v^n 5 ) (' ëïïwO’ 

_ . ? 0 — SïrzO’ 

A 7* 

et parce qu’on a déjà trouvé —,■=?,—, il en résulte le rapport cherché 
a' — i V 6n’Z’ 6n ’zV’ 

Observons de- plus que des valeurs précédentes, on déduit 

r 4- r° sa anz + TiyTrj ... 

r'+ r = alTZ'-f- n>'*Z'*; 

fl f't 

donc dans les termes , g - , on peut substituer les va- 
leurs suffisamment approchées tIZ=j(r+/'), nZ’=;(r+r'), 
ou simplement nZ=u , nZ'=u', en appelant u et u' les rayons 
vecteurs moyens et -ï(r'4-/‘), des arcs C’C et CG' 

parcourus dans les tems 8 et 8', on aura donc enfin 

a __ 9 / »' , V'\ 

7*"H' 6» i + 67v» 

rapport qui doit être exact aux quantités près du quatrième ordre ; 
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et parce qu’on a suppose A-f-A'^ i, il en résulte 

**r£»D+«r» (£-£>)]• 


3. Pour avtfir maintenant la valeur de e qui fait connaître 
celle de SD=#r(i — e), j’observe que les triangles C'SD, SC°C', 

donnent SD= AÀ - . , sin C“ = jr sin (S + £')■, donc 


SD = A/. Mais on a = ^+p} L~JrÙ; 

Mu b imt y(i -f- y *) 

A/' = -^i+p— } donc SD = et 

dCy-f-y ) O -hr/Wy i+yy«* v J 7 

par coasequcat 

e= _ yy / ('+ ^) 

‘ +w'u* 


Si l’on néglige les quantités du quatrième ordre, on aura sim- 
plement c=j.y(i+u*), ou en substituant la valeur de jÿ, 

c — n - Dans cette formule, où l’on doit rejeter 

les quantités du quatrième ordre, il suffira de faire nZ= [{o-f-r 0 }, 
ou même nZc=|/(rr*), car la différence c-H' — a \/ (rr°) est une 
quantité du second ordre; on fera semblablement 11// = V^( r ^)» 
et on aura enfin sous une forme très simple 


89 ' • 

résultat qui doit être exact aux quantités du quatrième ordre 
près. 

4. Les formules précédentes offrent un moyen facile de dé- 
montrer le théorème de Lambert, dont nous ferons usage dans 
cette nouvelle solution. Considérons, pour cet effet, le secteur 
C’SC décrit dans le tems fl; soit la corde CC ° = c°, puis 
r- (- r* -f- c* = a/> , r-f-r" — c°= 2 q,cc qui donne r-{-i‘s=p-\-q t 
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f ’=p — (j. Par le triangle SClfl, on a c**s=r'+r M — arr*. 

; =(''+0*-(7TO ; donc (TTV^=^ 0U ^ 7 ) = ^ 

= n .i-±ü‘ = nZ; donc z = ^. 

I •+■ yu ' n 

D’un autre côte, on a r* = ^ 1 y .^ , r+r* ou*>-f-<7 = — ■2^— 

’ (l+VU) T" ■ / J^-yU. 

Dy»/ | J. U *\» % 

+ ~ (V-f-Tû)‘ ~ > donc / , + 7 = aV / C/ , 7) + n 5'‘Z.*, « qui donne 
y Z= ,y ~ ; substituant cette valeur , ainsi que celle de Z , dans 

T 11 • 

l’équation du tems = yZ‘ -f- i> 3 Z 3 , on en tire 

8 = y O*— 7*) = l (** + r-f- c*) T — £(»* -f - r — c ’)*i 


c'est le théorème de Lambert qui donne l’expression du temps 
parla somme des deux rayons vecteurs et la corde de l’arc décrit. 

5. Ces différens résultats supposent que l’orbite est parabo- 
lique*, cherchons maintenant les valeurs des quantités A, A', e, 
lorsque l’orbite est elliptique ; ces valeurs deviendront applicables 
au secteur décrit par la Terre dans le temps 0 + 0'; et en suppo- 
sant le grand axe infini , elles devront coïncider avec les résultats 
précédens. — * — - 

Désignons par x,y, s, r, les coordonnées et la distance 
du point C , rapportées au centre des forces S, pour une époque 
déterminée t=ji désignons semblablement par x', y, a', r J t 
celles du point C' pour une époqne peu differente t=j- f— ô', 
et par x’, y, z', r*, celles du point C*, pour l’époque t=j — 0. 
Il s’agit de développer ces dernières quantités, suivant les puis- 
sances de fl' et 6, en supposant que les intervalles de temps fl' 
et 0, mesurés par les arcs correspoudans du moyen mouvement 
de la Terre, sont assez petits pour qu’on puisse négliger les termes 
du cinquième ordre et au-dessus. 

Ce développement se fera, comme dans l’art. 5 du premier 
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Mémoire, d’après les équations différentielles du mouvement; 
il en résulte, pour les coordonnées x' et x% des expressions de 
cette forme 

*•= x A + J «•(-«)» 

dans lesquelles les fonctions f el <p se développent jusqu’à la 
quatrième puissance de r inclusivement, comme il suit: 


A r ) = 1 — 53 + ^ — t 4 (gjr — g? + g^s), 

ç(r)=r -£ + £ 


Les cocfficiens de ces fonctions dépendent, comme on voit, 
du rayon vecteur r, de la quantité k = relatifs l'un et l’autre 

à l'cpoque <=y, et du demi-grand axe a de l’orbite supposée 
elliptique, lequel devient négatif dans le cas de l’orbite hyper- 
bolique, et infini dans celui d’une orbite parabolique. 

6. Les memes formules s'appliquant semblablement aux deux 
autres coordonnées, on déterminera, d’après les intervalles de 
temps 6' et 0, la position des points C' el C°, comme il suit: 

r=jr/-®)+3r*M), 

* = r/(-0)+ g *(-«) , 2' = 


Si l’on élimine de ces équations les coefficicns différentiels^, 

fS. — . on en déduira les trois suivantes: 
dt ’ <* ’ 

• ar*9(fl') — jt'^( — S) = anv , 

— j>{— 6) =^iv , 

• ■ z^8')— s'<p(--«)=2if, 

3 . 
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où l’on * fait w=aj{ — 0)|>(8') -*-,/$')$(■*■*■ 8)* c’est-*-difc en dé- 
veloppant jusqu'aux quantités du quatrième ordre inclusivement, 

w = a + r - + p (8+0') s (0' - 0)- 

On pourrait, des trois équations précédentes, éliminer entièrement 
p(0'), p( — 8), et tv, ce qui donnerait 

—/*) +j(~'x' — * V) ■+• z(x-y — x'y) = o i 

c’est la condition qui exprime que les trois points C% C, C', 
sont situés dans un même plan avec le point S, origine des 
coordonnées. 

7. Maintenant, de ce que les deux parties C°D, DC', de la 
corde C’C', sont entre elles comme A est ù A', il suit que 
leurs projections Q°E, EQ', sur le plan de l’écliptique, sont 
également entre elles— it A: A'; et comme les coordonnées des 
points Q% E, Q', dans le sens des x, sont x°, (1 — e)x } x', ou 
■nra 

A[x' — (1 — e)x]= A'[(i — e)x — jc 0 ], 

ou 

Ax' + AV = (1 — e)x; 

donc 

O 'A- 6 ) -t- MA'» 4- [A'<p(- 8 ) + AP(S')] £ =(, —e)xj 

on aurait semblablement, dans le sens des jr t l’équation 

[AX- 0 ) + A/ 0 ')] 7 4- [A'<P(- 0 ) + Ap( 0 ')] |=(1- e)j. 

Ges deux équations, où les quantités , sont 

dantes de 6 et ô', ne peuvent subsister à moins qu’on n’ait 

1 _e = A'/- 0 ) + A/{&'), 
o = Ap(— 8) -f- Ap(6'). 

La dernière fait connaître le rapport de A à A', lequel est 
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A' 


—çrr 3 * 


9* W 
9 

{F fcS^î 

I ’"ÏP + 5* 


et puisqu’on a A = ^, et r* = r— -fi '^= r— ^ il s'ensuit 

que si l’on néglige les quantités du quatrième ordre seulement, 

#* £9> 9‘ > . i 

on pourra faire i— ^ — ^=1 £- -» : on aura de meme 

6(r*r)* 

i — £“1 + 73 = I 1 î donc la valeur de -7 se rédqit à 

br 6frO* 

cette forme très simple . 

1 * 

1 — ■ -« 

A _ I G (r»r)* . 

a' ^ r • ' — ~/ï~ » 

1 1 

6(n') i 

et puisque A-j-A'= 1, on en déduit 



*ffT 

6(r»r)* / 


8. Substituant ces valeurs ainsi que celles de _/{0’) et de /"(“S), 
dans l’équation 1 — e = Xf{ — 0)+A/£fi’) , on aura, en négligeant 
toujours les quantités du quatrième ordre. 


ü A89'(9'— S) fly 
aï 3 ar 5 a rW 


Ainsi on voit que les valeurs de A, A', e, qui ont lieu dans 
l’hypotbèse d’une orbite elliptique ou hyperbolique , s’accordent 
avec celles que nous avons trouvées directement pour l’orbite 
parabolique, attendu que le demi-grand axe a n’entre point dans 
leur expression, au moins lorsque l'approximation n’est poussée 
dans ces valeurs que jusqu’aux quantités du troisième ordre ia- 


ds 
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clusivement. Pour avoir un degré d'approximation de plus, il 
faudrait que ]e développement des fonctions f( t), ç( t), fut 
continué jusqu’aux termes du cinquième ordre , ce qui , au reste, 
est superflu pour la solution de notre problème. 

g. Dans les applications ordiuaires, on pourra se borner aux 
termes du second ordre , ce qui donnera 

, t , tï — _ _ H' 

* + 7+7 \ 6H )*' e ~~ ar” 



• ai* f'* 

mais nous ferons remarquer que la quantité 1 — j aya nt 

, . . , G(rr»r 6( n'y 

cte mise dans les valeurs de A et A pour 

T?~ 4? ’ 

— #* 

on ne peut la réduire au premier terme > qu’en com- 
mettant une erreur qui sera quelquefois aussi grande ou même 
plus grande que le terme conservé. 

Le même inconvénient n'a pas lieu sur la valeur de e; car la 

quantité étant mise pour — — — ^F~ le second terme 

sera toujours beaucoup plus petit que le premier. 

Il est inutile d’ajouter que dans les premüra* approximations 
on pourra, en négligeant les quantités de second ordre, faire 
simplement 

r+v> = ¥+j ' > e = °- 

10. Venons maintenant aux équations nécessaires pour la so- 
lution du problème. La relation entre les coordonnées des points 
C°, C, C', fournit, ainsi qu’on l’a déjà vu, les trois équations 
Xx 1 +• X'x“ — (t — e)x = o, 

CO . — (i — e)jr == o, 

A s! -f--A'a* — (i — e)* = o. 
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Mais, par les donne’es de l'observation , on a, pour le point C, 

x — <r cos a + R cos A = <r cos a -f- X, 

( a ) J =■ * s *n a R siu A = a sin a 4- Y, 
a == a tang b ; 

et comme les points C*, C', fournissent de semblables valeurs, 
on en dc'duira les trois équations 

Au' cos a' -f- AV cos a" — f i — e) a cos a -f- w cos A r= o, 
A^ sin «' -f- AV sin rf* — (i — e) a sin a asr sin At= o, 
A* 7 tang A' -f- AV tang A* — (i — e) g tang b = o, 
où l’on a fait pour abréger 


AX' -f- A'X° — (i — e)X = «r cos A i , 

AY' -f- A'Y' — fi — e)Y = tr sin Ai. 

✓ 

Ces deux dernières quantités étant relatives au mouvement de 
la Terre, on pourra les réduire à une forme plus simple, au 
moyen de l’analyse précédente; car, en conservant l’analogie 
des dénominations, on a 


x-=XF(_e) + ^*(-fl), 

X' = XFCÔ')+-f 


F(-ô)=. 

Fceo =; 


JS* K fl 1 

aR 3 5R3> 
6 > KS’ 3 
aR 3 aR 4 > 


®(- 0 ) = -0 + ê^,+ 


Kf* 

4R 5 » 

i'— 6 '' _i_ ü 6 ! 1 
6R’ 4R5' » 


de là résulte, en ne négligeant que les quantités du quatrième 
ordre, AX' -t- A'X° — (i — e)X, ou 


_* vV * , , dX W, n , , N 

«A.=Xfe 5 r5+ —Sût - ) + ~ t ■ 


Soit encore par analogie 

t, 84' 


K88- 


sr 
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on aura 

XX' + A'X« - ( i - e)X = (e - E) (x + tfi . . 

enfin, si on appelle Xi ou Ri cos Ai, la valeur de X corres- 
pondante au temps y+i(8' — 6)» qui lient le milieu entre les 
temps des trois observations y — ô, y, y 4 - fl', on aura pour der- 
nière réduction 

XX' + X'X a — ( i — e)X = (c — E)X i = (e — E) R i cos A i . 

On trouverait semblablement 

XY'+X'Y* — (i — e)Y=(e — E)Y i = (e— E)Ri sin A i ; 

nos trois équations se réduisent donc aux trois suivantes 

o=X(7 J cosn'+X'fl'*cosa* — ( i — e)rcosa-\-(e — E)R i cos A i , 

(3) osa^.tr'sino' -f-AVain o° — (i — e)<rsin «-(-(e — E)Risin Ai , 
o=X9 / tangi'-j-AVtang b " — ( i — e)irtangê. 

C’est de ces équations que dépend principalement la solution 
du problème ; mais elles demandent à être employées avec 
quelques précautions, pour qu’on puisse en tirer des conclusions 
certaines; il est donc nécessaire de faire, sur ces équations, 
plusieurs observations préliminaires qui serviront à diriger les 
calculs. 

ii. La première observation que présentent les équations 
(3), c’est que deux des inconnues e - *, y, a 1 , ne peuvent pas 
en être éliminées à la fois. 

En efTet, si on élimine et a' de ces trois équations, on en 
déduit 


. e — E ^ tin (A i — o°)tanR6'— .»in(Ai — a' ) tangi* 

(‘0 V , — e 1 " jin(a — a°)taugi' -f- ûn(a' — a) langé" -J-sin (u° — u'Jtangl’ 


or, qnand même on connaîtrait d’avance la valeur de e, qui 
suppose r connu, on ne peut pas regarder comme suffisamment 
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approchée la valeur de a donnée par cette formule , attendu 
que le dénominateur est, comme nous l'avons déjà fait voir, 
une quantité du troisième ordre, susceptible d'être considéra- 
blement affectée de l’erreur des observations. 

Nous donnerons, par cette raison, le nom d 'équation douteuse 
h l’équation (4), et nous nous abstiendrons en général d’en faire 
usage, pour déterminer la valeur de a. Le même nom con- 
viendrait aux deux équations semblables qui serviraient à déter- 
miner et cr'. 

Comme la supposition de l’orbite parabolique donne une équa- 
tion de plus que d’incoonues , on pourra toujours parvenir à 
la solution, sans employer l’équation douteuse. Il y a néanmoins 
quelques conséquences à en tirer. 

12 . Si le dénominateur de la formule (4) est nul, c’est-à-dire, 
si l’on a 

D=sin(« — rt°)tangi'-(-sin(<j' — a)tang£ a -f-sin(a* — «')tangé=o, 

auquel cas les trois lieux apparens de la comète sont situés 
dans un même grand cercle, il faudra que le numérateur de 
celte formule se réduise aussi à zéro, ce qui peut avoir lieu 
dans deux ans. 

i*. Si l’on a 

P = sin(Ai — d°) tang é' — sin (Ai — «')tang£’=o, 

c’est-à-dire si le lieu du Soleil au temps /-f-f (V — 8), moyen 
entre les temps des trois observations, est situé dans le même 
grand cercle que les trois lieux apparens de la comète , alors la 
valenr de cr devient J, et il n’y a rien à en conclure. 

a*. Si la quantité P n’est pas zéro, ou si le lieu du Soleil 
n’est pas dans le grand cercle qui passe par les trois lieux ap- 
parens de la comète, mais qu’il soit très près de l’une des in- 
tersections de ce cercle avec l’écliptique, il n’y aura encore rien 
à conclure de l’équation (4). 
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Mais si la quantité P, qui est généralement du premier ordre y 
est assez grande pour que son signe ne puisse pas être changé 
par les erreurs des observations, alors il faudra qu'on ait à très 

peu près e — E = o, ou p = , ce qui apprend que la comète 

et la Terre sont également éloignées du Soleil au temps de l’ob- 
servation moyenne. 

Cette circonstance facilitera les premiers calculs d’approxima- 
tion; car, sachant qu’on a à très peu près r=R,on devra avoir 
aussi à peu près 

«■ = aR cos* b cos (180* -f- À — a) , 

180* + A étant la longitude du Soleil à l’instant de la seconde 
observation. 

Telle est donc la conséquence qu’on pourra tirer de l’équa- 
tion douteuse , lorsque la quantité D étant nulle ou très petite , 
la quantité P est en même temps assez grande pour que son signe 
ne puisse pas être changé par les erreurs des observations. 

îa. En général, lorsque les quantités D et P sont assez grandes 
l’une et l'autre pour que leurs signes soient bien déterminés, 
ce qui est la circonstance la plus favorable à l’exactitude de 1 équa- 
tion douteuse, on conclura de cette équation que le signe de 

est le même que celui de e — E, ou de p — on saura par 

conséquent si r est plus grand ou plus petit que R. 

11 pourra même se présenter des cas où l’exactitude de l'équa- 
tion douteuse sera suffisante pour qu’en combinant celte équation 
avec celle que donne le triangle SCT, on en déduise des pre- 
mières valeurs approchées de a et r. 

En effet, si l’on fait , | = h, l’équation ( 4 ) pourra être 

mise sous la forme 

7 = h + hr '> 


Digitized by Google 


( «7 ) 

d’ailleurs le triangle SCT donne 

r* = R* + aRff cos (A — fl) -H 


C05*t* 


Ces deux équations sont satisfaites, lorsque a=o et ;-=R; 
mais elles auront une autre solution, si cos (A — fl) + j/dl* est 
une quantité négative, et celle solution s'obtiendra très aisément 
après quelques essais. 

Le cas que nous venons d'indiquer ne pouvant être que fort 
rare, nous devons en faire abstraction, et regarder en général 
les équations (3) comme ne pouvant servir à éliminer à la fois 
deux des inconnues a°, or, a'. Il reste donc à examiner seule- 
ment les combinaisons dans lesquelles il n’y a qu’une inconnue 
d'éliminée, et qui serviront par conséquent à déterminer deux 
des inconnues e*, a - , <f, jpar la troisième à laquelle on donnera, 
dans les premiers essais, une valeur hypothétique. Ces combi- 
naisons peuvent être présentées sous beaucoup de formes di- 
verses; nous allons rapporter les principales, afin qu’on soit en 
étal de choisir dans chaque cas particulier les formules les plus 
propres à donner des résultats exacts. 

j 3. La première combinaison qui se présente est celle des 
deux équations suivantes: - J — - — - — ' 

♦ 

X»'sin(Ai — fl')-f-A'ff*sin (Ai — fl*) — (i — e)ffsin(Ai — a)=o, 
Ao 7 lang b' + A'ff” tang b‘ — (i — e)a tang b : 

iùa 

d'où résultent les formules 


:OJ 

# ÿ*i*. y 


( 5 ) 


^ (î — (Qj- ^taiigê'ai n^Vi — a ) — tangfc sin(A i — a) 

V tangi'sin (Ai — a") — t*ngè“sm(A I — a')’ 

, __ (i — e)r tangi »in(Ai — n* ) — -tang&»sin(Ai —a ) 
“ A ‘ tang6'sin(Ai — u°) — tang4“sin (Ai — a) 


On pourra , dans les premiers essais , faire e = o, A = —-y, , 

rT® 

A' = > et alors ces formules donneront aTec une exactitude 

«+» » 

3 
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suffisante les rapports — , —, pourvu que leur dénominateur déjà 

désigné par P, ne soit pas une quantité trop petite. Or nous sa- 
vons que ce dénominateur est nul, lorsque les lieux apparens 
de la comcte, dans la première et la troisième observation, sont 
situés dans un même grand cercle avec le lieu du Soleil à l’ins- 
tant de l’observation moyenne, ou pour mieux dire, à l’instant 
qui lient le milieu entre les temps des trois observations. 

Toutes les fois donc que celte coïncidence aura lieu ou exac- 
tement ou à fort peu près, il faudra recourir à d’autres for- 
mules , mais le plus grand nombre des cas pourra se résoudre 
par les formules précédentes dont nous avons indiqué l’usage 
pour la première approximation. Les approximations suivantes 
exigent qu’on substitue successivement des valeurs plus appro- 
chées de A, A', e , ainsi qu’il sera expliqué ci-après. 

14. Une seconde Combinaison qui produit un résultat très 
simple , consiste à n’employer que les deux premières des équa- 
tions (3); on en tire les valeurs de a‘ et <r', exprimées par le 
moyen de s - , comme il suit : 


(G) 


(1 — e)r $in ( a ' — a) -f- (E — e) R l sin (a' — A j ) 
a' sin (a' — a°) 

(1 — e)c-siti ( a — cr) -f- ( e — E)R i «in (a° — Ai) 
a sin (a* — a 0 ) 


Dans les premières approximations, on pourra supposer E — e=o, 
1 — e— 1 , A = , A' = —y, ce qui donnera simplement 


1 + 1’ rsin ( a — à) S — j— S' rata (rc — o») 

7 " si n (a' — a”) ’ « sin (a — a*)’ 

Après les premiers essais, il faudra substituer les valeurs plus 
exactes de A, A', e, qui ont été données ci-dessus art. 9. 

Les équations (6) s’appliqueront avec succès toutes les fois 
que le mouvement en longitude sera assez grand par lui-méme 
et plus grand que le mouvement eu latitude. Les équations ( 5 ) 
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seraient préférables, dans le cas contraire, sauf l'exception dont 
nous avons déjà parle'. 

i5. Pour donner un exemple de quelques autres combinaisons 
qui peuvent offrir des avantages particuliers, multiplions les 
équations (3) respectivement parcosip, sin p, — cotè’cos(p — a'), 
nous aurons, en ajoutant les produits, celle nouvelle formule 


(i — — a) — tan£icr>tiVos(0— o 0 )]-f-(E— e)Ri®ns(£ — Ai) 

cos(^— a')— langt'cotï°coi(0 — a°) 


dans laquelle <p peut être pris à volonté. 

Si l’on fait successivement <p = o, £= 90 *, on aura les deux 
formules qui suivent, pour déterminer a' par le moyen de v 


(7) 


(1 — f)r(cm a — Iang6cnt6°coi 0°) -J- (E — e)Ri COJ Al 
A (cos o' — lang A' cot A' cos a") 

(1 t~ s) r(sin a — tan pAcnt A*»ina*) -f- (E — e)Ri sin Al 
A (siu a — lang A' cul A 4 * aina*) 


On choisira, pour chaque cas particulier, celle des deux for- 
mules qui aura le plus grand dénominateur. 

Les formules, pour déterminer a’ par le moyen de cr, se dé- 
duiront des deux précédentes, par la simple permutation des 
signes ' et •; on obtiendra ainsi les deux formules suivantes 

(1 — f)r(c ns a — tang A cot A ' cos ai ) -f- (E — <*) R i cos A 1 

,o\ a' (cos o a — taug If cot b' cos a') * 

' ' ( 1 — e) r (ain o — lang A cot b' sin n') - 4 - (E — «^RisinAi 

a' (»in a° — tang A" cot A ; sin «') 

El on voit que ces formules ne sont guère plus composées que 
les formules ( 6 ). 

16 . Ayant montré l’usage qu'on doit faire des équations (3), 
suivant les diflerens cas, nous allons rassembler sous un même 
point de vue les diverses opérations nécessaires pour résoudre 
les équations du problème, et eu déduire les élémens de l’orbite 
parabolique avec toute la précision que les observations com- 
portent. 


1 
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I. Les premières données du problème sont les trois longi- 
tudes de la comète a°, a , a', et ses trois latitudes b’, b t b‘, 
observées à des époques j — 0, j , dont les intervalles 

sont ô et 0'. Ces quantités 0 et 0', données en jours moyens, 
devront être multipliées par le nombre qui a pour logarithme 
8.23558ai , afin d’être converties en arcs du mouvemeat moyen 
de la Terre, ces arcs ayant pour rayon l’unité. 

Pour plus d’exactitude, nous supposons toujours que les lieux 
donnés de la comète, sont des lieux observés, et non des lieux 
déduits de l’interpolation. 11 n’est pas nécessaire que les inter- 
valles de temps 0 et 0’ soient égaux, mais il est bon qu'il n’y ait 
pas entre eux une trop grande inégalité, et que le plus grand 
des deux n’exccde guère io jours, parce que la solution ana- 
lytique est fondée sur ce que les quantités du quatrième ordre 
en 0 et 8', peuvent être négligées. 

On calculera par les tables astronomiques, ou bien on prendra 
dans les éphémérides les lieux du Soleil pour les iustans des 
trois observations. Ces lieux augmentés de 180° donneront les 
longitudes de la Terre, A% A, A'. On calculera de même les 
rayons vecteurs eorrespondans R*, R, R'; il faudra ensuite en 
déduire les valeurs de A et de R pour le temps r (ô' — 0) 
qui lient le milieu entre les époques desTïotSTobservations. Ces 
valeurs, désignées ci-dessus par Ai et Ri, se trouveront avec 
une exactitude suffisante pour l’usage qu’on en doit faire, par 
les formules Ai =^(A“-1- A-f- A'), logRi = ^log(R”RR'). On 
pourra cependant , dans le dernier calcul d’approximation, sub- 
stituer aux quantités (e — E)RicosAi, et (e — E)RisinAi, les 
valeurs de •w cos Ai et ■arsinAi, dounées dans l’art, io, les- 
quelles sont indépendantes de toute supposition sur la figure de 
l’orbite terrestre. Mais cette correction qui ne laisse pas d’allon- 
ger le calcul, est en général de peu d’importance, parce quelle 
est de l’ordre des quantités qui ont été négligées. 
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II. D’après les données qui viennent d’être mentionnées, on 
pourra calculer les valeurs des quantités 
P=langê'sin(Ai — a*) — langi*sin(Ai — a'), 
D=tangi'sin(a— a°)-(-tangê°sin(n'-— a)4-tang£sin(a* — a'); 

si elles sont assez grandes l'une et l’autre pour que leurs signes 
ne puissent pas être changés par les erreurs des observations , 
on en conclura que R — r est du même signe que DP; et si en 
même temps cos(i8o° + A — a) est positif, on voit par l’équation 

r* = R* + aRor COS (À — a) -f- —, 

que aR cos (180"-}- A — a) — s devra avoir le même signe que 
R — r, ce qui fournira une première donnée sur la valeur de 
Si D est très petit par rapport i P, on en conclura que les 
distances r et R sont à très peu près égales. 

Si D et P sont très petits l’un et l’autre, il n’y a plus rien 
à conclure sur la grandeur relative de r et R , et on tombe alors 
dans les cas d’exception dont nous avons parlé. 

Ifl. Cela posé, on fera une première hypothèse sur la valeur 
de tr, soit arbitrairement, soit d’après les résultats du calcul 
précédent; ensuite on calculera les valeurs de a° et a', par les 
formules (5) ou par les formules (6), ou même par les formules 
(7) et (8), suivant les différons cas. 

Dans les premiers essais, on pourra, en négligeant les quan- 
tités du second ordre, faire c—o, E=o, A= ■ ■ ■ * ■■■■ ; , A' = , 

ce qui simplifiera d’autant les formules; on parviendra ainsi à 
des premières valeurs approchées de a, a'. 

Celle première approximation* étant obtenue, il faudra, pour 
aller plus loin , substituer les valeurs 
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alors les formules qui expriment a • et a contiendront à la 
fois ir et r, mais la marché du calcul n'en sera point arrêtée; 
car ayant fait une première hypothèse sur la valeur de <r, on 
en déduira immédiatement r par l'équation que donne le triangle 
SCT, savoir 

/•* = R* 4- 2R.tr cos (A — a) + ~j, 


et substituant ces valeurs dans les expressions de ff* et </ , on 
connaîtra ces deux quantités. 

IV. On calculera ensuite les deux rayons vecteurs correspondans 
r* et r par les formules 


(10) 


*" = R** + îRV cos (A" — «•) + — .-p. 


r- = R'* + 2 RV cos (A' — a') + ~rp- 


Quant à la corde c qui joint les points C", G', on pourrait la 
calculer par la formule 


. r . /* 
e* = — r 

A ‘ A 


(l— 0* 

A\' 


mais, comme cette formule suppose qu’on connaît exactement 
les valeurs de A, A', e, ce qui ne peut avoir lieu que lorsque 
l’approximation est déjà fort avancée, il conviendra de calculer 
c par une autre formule qui est plus composée, mais qui n’est 
sujette à aucun inconvénient. 

Cette formule n’est autre chose que l’équatiou 

C * = (>•' — x°y + (/ —yy + (*' — 

dans laquelle il faudra substituer les valeurs 


x? = R° cos A 0 -f* 9* cos a° t y 9 = R°sin A° -f- r° sin a° , t° = r 1 tang b 9 , 
x = K cos A'-jV cos a', y' = R' sin A' -f- r sin a, z ~r lang b', 

ce qui donnera, en introduisant une auxiliaire V, 
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{ c , =r« + r / *— aV, 

V = R°R' cos (A' — A 0 ) -f- RV cos ( A' — o') -f- «®c coj (a' — o*) 
+RV cos (A* — a) -f- «V tang b° tangi' ; 


enfin pour ve'rifier toutes ces valeurs, on a l'équation des temps 
donnée par le théorème de Lambert, savoir: 

(.a) 6(8 + G') = (/* +/' + *)* — (/- + ✓ — c)‘; 

« 

c’est au moyen de cette dernière équation qu’on rectifiera suc- 
cessivement les différentes hypothèses faites sur la valeur de a, 
et qu’on approchera de plus en plus de la vraie solution. 

V. Lorsqu’on sera parvenu à une valeur de or exacte dans les 
deux ou trois premières figures, on connaîtra en même temps, 
avec un semblable degré d'approximation , les valeurs des rayons 
vecteurs /*, r, r 1 . On sera alors en état d’évaluer avec un degré 
d’exactitude suffisant pour la dernière approximation, les quan- 
tités e. A, A', ce qui se fera par les formules 




m<'* 


6(rV)* 6(r/)^ * 


(. 3 ) log a — i°g (jr^rf) ~ r+7 C» 
log a' = io g -i- C- 




On connaît déjà la quantité E = R * R - j p ; ainsi en substituant 


ces données dans les équations du problème, on parviendra à 
satisfaire à ces équations aussi exactement que le comportent 
les observations soumises au calcul. 

VI. Les équations étant ainsi résolues, il ne reste plus qn’à 
déterminer les élémens de l’orbite. C’est ce qu’on fera aussi 
simplement qu’il est possible par les formules suivantes. 

Etant donnés les deux rayons vecteurs extrêmes r”, r', et la 
corde c qui joint leurs extrémités, il faut d’abord trouver l'ano- 
malie 4‘ *l u * répond au premier rayon vecteur r’ et la distance 



( 2 4 ) 

périhélie fl; soit pour cet effet s=-j(r*-f-r'-f-c), u=4(r*-f-/-'— -c), 
on aura les formules suivantes, où l’on suppose qu'à l’instant de 
la première observation, il reste à la comète l’anomalie 4° 3 
parcourir, pour arriver à sou périhélie. 


(»4) 

n 


sin f4* 

taug ’ 4° 


= r* cos‘j4° 



_ r°-\ /( su) 

(S-T°)(S-/) 

— Ws—ÿ U y • 



il sera bon de calculer 4* ainsi que fl, chacun par deux for- 
mules différentes, afin d’être sur du résultat jusque dans la der- 
nière décimale des logarithmes. 

Counaissant n et l’anomalie 4% on trouvera par la table des 
comètes K temps employé à parcourir l’anomalie 4’» ce temps 
s’ajoutera à l’époque de la première observation, ou s’en retran- 
chera, selon que 4’ sera positif ou négatif, et l’on aura l’ins- 
tant du passage au périhélie. Il sera même encore plus simple 

de calculer directement le temps par la formule ; 

1 _ 

1 = ( 2 n s )* (langi4-H lang 3 H) > a y an t soin d’ajouter au loga- 
rithme de chacune des parties le logarithme constant >, 7644 ‘ 79 * 
qui est le complément de 8,a355Sai , afin que t soit exprimé 
en jours moyens. 

Désignons par a* et S‘ la longitude et la latitude hélioccntrique 
du pointC*, par «' et €' celles du point C'; on aura pour dé- 
terminer ces quantités les formules suivantes où tout est connu, 
puisque les distances accourcies u° et <r' sont données par la réso- 
lution des équations du problème. 


(t5) 


sin£*=^ langé’, sin£ =^ tangê'. 


sin(a«-A*)=[-^sin(a°-À*), sin(*'-A')=|^sin(«'-A') 


tanf^ 
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Le mouvement sera direct ou rétrograde, selon que a.'—**, 
sera positif ou négatif. 

Soit maintenant N la longitude du nœud ascendant, et I l'in- 
clinaison de l'orbite, on aura la double équation, 

(16) cot I s= col sin («• — N) = col C sin («' — N), 

d’où résulte d'abord pour déterminer N la formule 

(> 7 ) <ang — N 

Dans l’usage de ces formules, les latitudes boréales sont tou- 
jours supposées positives et les latitudes australes négatives; 
on suppose aussi l’inclinaison 1 constamment moindre que 90°, 
de sorte que cotl doit toujours être positive. Il faudra doue 
que sin (*•’ — N) soit de même signe que la latitude c’est 
ce qui fixera le choix entre les deux valeurs de IV qui corres- 
pondent à une même valeur de tang ( * — ■ N). 

Enfin si l’on appelle II l’angle compris dans le plan de l’orbite 
entre les deux rayons vecteurs menés au nœud et au point C% 

on aura ■ * - . ■ . « 

(18) w cos H = cos 6* cos (a* — N); 

et si à la longitude du nœud on ajoute la distance du périhélie 
au nœud 4 * + H, on aura le lieu du périhélie sur l’orbite 

= N + 4.° + H.* 

Application de la seconde mclliode a la première comète de 1 8 o 5 . 

16. Nous choisirons dans le tableau de la pag. 14, premier 
supplément, les trois observations suivantes qui serviront à 
déterminer l’orbite. 


sin -f£2 
" a in (C' — O 1 ) 


tang 


( a6) 




Époques. 

Longitude*. 

Lat. bor. 

Lieu du 0 -i 8 o°. 

Log R. 

Octob. 33.68488 
30.68670 
34.73088 

i 63 ° 3 o' 53 " 
i 83 . 48 . 3 a 
191 .46. i 5 

aa’ 5 g' 53 " 
15.37. 3 ) 
la. a.ag 

aq°io' 5 o“ 

37.19.48 

4 1 . aa . a6 

9-997431 
9 • 9964 f )6 
9.996053 


Dans cet exemple, on a ê = 8'ooi8a, fl' = 4'o34 l8 »* . . 
6 — j— Q = i a'oTiGoo ; ces nombres étant multipliés par un factcnr 
dont le logarithme est 8.a3558ai, on a les résultats suivans , 
exprimés en logarithmes 

0... 9.1387708 0'... 8.84i3373 

0 + 0 '.. 9.3160643 0 + 8 '... 9.3160643 * 

r+7-” 9- 8aa 7° 65 r+V” 9- ^aySo. 

La seconde colonne du petit tableau qui précède, contient 
les valeurs des quantités a°, a, a!-, la troisième, celles des quan- 
tités 4°, b, b'; la quatrième, celles des quantités A", A, A'; 
enfin la cinquième, celles des logarithmes des distances R°, R, 
R'. Si l'on prend le milieu entre les nombres de ces deux der- 
nières colonnes, on aura, avec une exactitude suffisante, les 
valeurs des quantités que nous avons appelées Ai,Ri, savoir: 

Ai = 56”o'4i", 
log Ri == 9.9966567. 

Pour le calcul des quantités P et D, nous rassemblons ici plu- 
sieurs logarithmes constans qui seront utiles daus les opérations 
ultérieures 


tar.g b ° 3.6378110 
langé 9.446581s 
tang b ' 9.3090147 


$in (u' — a) 9.1414978 sin(a' — Ai) g. 6 i 33857 
sin («' — a”) 9 . 677583 1 

sin (a — a °) g. 54353 o 5 si 11 (a 3 — Al) 9.9004139. 


Par ces valeurs, ou voit immédiatement que la quantité P 
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est positive et assez petite; si on la calcule exactement, on 
trouve P = o.oo46585. Quant à la quantité D, si l'on convient 
de désigner par [N] le nombre dont le logarithme est N , on 
aura 

D = [8.87a545a] -f- [8. 7693388] — [9. 134184^] =0.0202617. 


Celte quantité est positive et fort petite; cependant son signe 
ne pourrait changer qu'en supposant des erreurs de 4 h 5' au 
moins dans quelqu’une des positions observées; et comme ces 
erreurs sont peu probables, il y a lieu de croire que les quan- 
tités P et D sont de même signe, et qu’ainsi on a la distance 
r<R, c’est-à-dire qua l’époque de la seconde observation, 
la comète devait être à une moindre distance du Soleil que la 
Terre. C’est toute la conséquence que nous tirerons de l’équa- 
tion ( 4 ); car d’ailleurs nous aurons occasion de faire voir que 
cette équation est fort défectueuse. 

17. La petitesse de la quantité P qui sert de dénominateur 
aux formules (5) , nous avertit que cet exemple tombe dans 
les cas d’exception; c’est pourquoi nous devons recourir aux 
formules ( 6 ) dont l’usage ne présente aucun inconvénient , 
puisque le mouvement en longitude est assez considérable. 

Ces formules sirapliGécs en faisant c = o, Eao, X=j-“, 
X' = 7 - 7 - 2 , se réduisent aux suivantes 

= <7 [o,o453^09], <7° = <7 [ 9 . 9586417 ], 


mais l’approximation qu’on en déduirait serait fort imparfaite , 
parce que, dans cet exemple, les intervalles 0 et fi' sont trop 
éloignés de l’égalité. 

Dans tous les cas seniblables , il convient de ne point se 
servir des équations ainsi simpliGées, et de substituer immé- 
diatement dans les formules ( 6 ), les valeurs de 0 et 0 ', tirées 
des formules (9), ce qui donnera les valeurs suivantes de a' 


( 38 ) 

el exprimées par les deux indéterminées a et r. 

" f = «£6. 93864 17] — p [7.740864] — ^[8.086264] +[8.0967763, 

<f = r[o. 0432409] — ~ [7.643860] + [8.075856] — [8.08G371]. 

On prendra d'abord une valeur hypothétique de *, d’après . 
laquelle on calculera le rayon vecteur r par la formule 

r* = 0.9839930 — a [0.2185267] +- e‘ [0.0336962]; 

r étant connu à peu près ( car les premiers calculs n’exigent 
pas beaucoup de précision ) , on calculera les quantités 0° et *' 
par les deux formules précédentes. 

Ensuite on calculera r°, r 1 , V et c, par les formules 

1 

/" = 0.9881934 — tr*[o.i 4 o 35 g 6 ] -+*“* [0.0719354], 

r'* = 0.9819876 — a' [0.2363369] +- o - '* [0.0193349], 

V = o . 963404 a — ar* [9 . 7 1 99493] + . 986770 1 ] 

— *'[9.9766975], 
c* = /**■+ r'* — 2 V ; 

enfin on vérifiera l'hypothèse, ou on en connaîtra l’erreur par 
l'équation des temps 

6(0 +- 9 ') = (r° + r' +• c)‘ — (1* ■+- r ' — c) 3 S= 1 . 2422687.' 

18. Comme il serait inutile de rapporter ici des essais de 
calcul trop éloignés du vrai résultat, nous ferons pour première 
hypothèse * = 0.6 = [9.7781512], ce qui donnera d'abord.,., 
r* = 0.3797624 =[9.5795120], r— [9.7897560]; on aura en- 
suite par l’application de nos formules 


J 
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*"—[3 - 66g 5087] , ?'=[g . 838 g 8 co] , 

r ,, =o. 6 oo 3 ago =[9.7783893] /•sso . ago66 i 3 =[g . 4634 17 a] 

^=0.7748089 =[9.8891948] o.6co 3g9Q [9.7317086]— r 7 * 

“ r'= ro. 0391476 ' r° , +r ,, z=o. 8909903 

r °+r' =1 .3.3.9 5T4 I or - aV— 0.7 Ô37670 

= 0.3 70 4366 1 1 c*=o. i37aa33=[9 . i374a 7 8] 

a=i. 6843930 b*= o. 9 164888 [9.56871 3g]=e 

b=o.g 435 ao 8 1 • “695887 

1 . 3439687 
erreur + a 73 a °° 

Soit pour seconde hypothèse e^o. 58 =[g. 7634280] , on aura 
les résultats suivans 


'*= o. 38739890= [9 . 588 1 554 ] , r=[9-794 c 777]. 

*-’=[9, 6553750], r'=[g .8343310] , 


r“*=o. 6047739 =[9.7815930] J / *=o. 3976389=69. 47 36897] 

,'==0.7776711 =[9.8907965] 0.6047739 [g. 7368448 ] =/ 

1^= 0. 5455639 r**-f-7 /, =o.goa4ia8 

^'=‘.3333340 ^ .8 P 

0= 0.3600 963 * c’=o. 1396694= [9. 1 138075] 

a=T 7 Ê 8333 o 3 b = [3.5564187]=* 
b=o . g 63 1 377 1 . o 3885 1 1 

1 .3433687 
erreur — 34176 


11 résulte de ces deux erreurs qu’en faisant n = = 0.11 1 jg T 

la quantité qgi a pour valeur A dans la première hypothèse et 
B dans la'seconde, aura pour vraie valeur A — /j(B — A). On 
trouvera ainsi la solution suivante: 


a ' ’=: 9 - 6569463] , =[9.8906184] , c— [9.5577858]» 

<7 =[9. 7 65 o 65 i], /• =[9.7935971], 

<r'=[9. 8359598], '•' = [9.7363757]. 


Dans les cas ordinaires, on pourrait s’arrêter à ce premier 
résultat, et on en tirerait des éléniens de l’orbite suffisamment 


(5o) 

approchés. Mais, dans le cas présent, où il y a une grande iné-. 
galité entre les intervalles de temps 8 et 8', les quantités A et 
n’ont pas été déterminées assez exactement, et on ne doit 
considérer le résultat précédent que comme une première ap- 
proximation qui a besoin d’être rectifiée. 

19. Au moyen des valeurs trouvées pour r*, r, r', il faudra 
calculer les valeurs de X, X', i — c, e — E, par les formules 


E= 


AOL 

R*RR' » 


e = i- M - ? 


Tné* 




G(r*r)‘ 6(rr)* 

log A= log ( 4 ^) iogx'=iog( r -p 7 ) + j± ? • £ , 


et on trouvera les résultats suivans : 

i 

( = [y . 564 5 1 07] , (e— E) = [8 . 1 35855 a] ,' 
x =[9.8219510], (1— e) = [9.9919393], 
X'=[9.5a68i 13 ]. 


Substituant ces valeurs dans les équations (6), on aura de nou- 
velles expressions de 7° et a' qui ne renferment plus que <r d’in- 
connue, savoir: 

< r ° = ff [9-9»9°4*®] — [8.54 i5o33], 
g ' = a[o. 0359557] -+- [8.5334i 17]; 

si l’on prenait pour première hypotbèse la valeur d!*b ■ qui résulte 
de la première approximation, savoir « = [9.7650651], on trou- 
verait que cette valeur est encore loin de satisfaire aux équa- 
tions du problème, et qu’il faut l’augmenter notablement, ce qui 
• prouve que la première approximation n’est pas suffisante. 

La raison en est que pour la première approximation, on a 
supposé les rayons vecteurs /•”, /•' égaux à r, dans les valeurs 
de A, A', e; cette supposition ne produit qu’une différence assez 
petite dans la valeur de e; mais dans celle de £ qui sert à 
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trouver A et A', elle en produit une fort grande; elle donne 
en effet £ = - = [7.6290203], quantité' presque double 


de la valeur plus approchée £ = [7.5643107]. 

20. Soit donc pour i re hypothèse ff=o. 6 o 45 =[g. 7813963], 
on en déduira 


r*=[g . 6799602] , r' =[9 83 9 37 c 5 ] , 

r^rso. 5973194 =[9.7763067] » / *^ro.29o5o8i=[9.463i58a] 

1^=0.7728645 =[g.888io33] o . 5 q7.3 1 g4 [9.7315791]=/ 

/=o. 538 g 88 o r° , -t-r / ’=o. 8878275 

>’®+K=i . 3 i i 85 a 5 J aV =o. 736461 3 

c — o. 36 a 4447 % 1-1 313 c*=o. i 3 i 366 a=[g. 1 184836] 

a=i .67^3072 h*= o. 9250798 e=[g . 55ga4 1 8] 

b=o. 9494078 1.2413714 

1.3423687 
err. — 8973 

Soit pour seconde hypothèse e- = o.6o55 = [9. 7821141], on 
aura les résultats suivans: 

r =[9. 68o73o3], r =[9. 8400529], 

r**=o. 597106a =[9.7760516] 1^=0. 2903 o 65 =[ 9 . 4627072] 

r*=o. 7727264 =[9.8880358] 0.5971062 [9.7313536]=/ 

r = 0.5387082 /”+/*=o. 88731 27 

r"+/= i. 3 1 14346 J ^ . aV =o. 7556 59a 

c= o. 3628408 , a ' 1 c ®° c*=o. i 3 i 6535 =[g. 1 194334] 

a=i .6742754 b’sro^sSSgoa , <=[9.5597162] 

b=o.g485g38 1.2425188 

1 , 2422687 
err. -f-aSoi 

D’apres ces deux erreurs, on trouve que le milieu à prendre 
entre les deux valeurs A et B d'une même quantité, données 
par les deux hypothèses, est B-]-(A-— B) (0.218). On obtient 
ainsi les valeurs logarithmiques suivantes des rayons vecteurs 
r, r> : 

r 0 = [9.8880427], r=[9. 7886828], r'= [9.7514028]. 



( 3 a ) 

Ces valeurs sonl assez différentes de celles qu’avait données la 
première approximation; c'est pourquoi il conviendra de pro- 
céder à une troisième opération dans laquelle nous donnerons 
aux équations (6) un dernier degré d’exactitude. 

ai. Voici d’abord les valeurs rectifiées de A, A', i — e, qui 
résultent des valeurs précédentes de r, r'; 

A=[9. 8219155], A'=[g. 5 a 684 ig]» > — £=[9.9918! 88]; 

substituant ces valeurs dans les formules de l’art. 10, qui de- 
viennent dans notre exemple 

«wcosA 1 =A[g . 87 1 353 o]+A '[9 .9378/, 19] — ( 1 — e)[g . 8969484] , 
wsin A 1 =A[g . 8 1 Ga347]— {— Av'[ç) . 68748 1 9} — ( i — e)[9 .7792587], 

on en déduit 

«ar cos A 1 = [8 . 03990 20] , lang A 1 = [9 . 87 1 67 1 1 ] , 

•arsin Ai=[7.9h573i |, Ai= 36 “ 3 g' 2 o", 

<ar=[8. 1555982]. 

On voit que l’angle Ai, ainsi déterminé, diffère sensiblement 
de celui que nous avions pris pour sa valeur approchée, savoir : 

Ai = 36 'o' 4 t"- 

Substituant ces valeurs dans les formules suivantes qui tiennent 
lieu des formules (6) ; 

(1 — e)*-sin {a — » n) — ar si n (a' — A 1 ) 

“ x'siu (u' — a“) * 

(1 — f) rsin (a — Q°) -f- »r sia (o° — Ai) 

** “ X sia (a — a") * 

on aura 

o ° = a [9.9288916] — [8 5552427] , . 

<r' = a [o.o 3585 o 7 ] [8.5401948]. 

À l’aide de ces deux équations, ou ne tardera pas à obtenir 
le dernier degré d’approximation dont nos formules sont sus- 
ceptibles. 
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aa. Soit 1°. a = [9.7790000]; on aura les résultats suivans ; 


*■*=[9 • 6763045] , r'=[g . 8373359] , 


7~=°. 59837755 =[9.7759753] /*=o.a9i43oi5=[9. 4645196] 

^=0.77354866 =[9.88848760] ^*=0.59837755 [9.7353598] 


r' =0 . 53983346 


r “+ r, =«-3i338aia a*=a. 1694653 

c=o.36a46583 , J ’ Q -, 
I b*=o. 9575854 


a_i .67584795 

b=o. 95091639 


1 .3431769 
1.5435687 


°-8897977 
aV=o. 7584163 

c*=o.i3i38i5 =[9.1185343] 

<=[9.5593671] 


err. — g 1 8 


Soit a*. = [9.779*000]; on aura 


r =[9-6763131], «-'=[9. 8374379], 


»**=o. 5983468 =[9-7769530] 

r*=o. 77353877 =[9.8884765] 


^=0.53979365 

f*+r , =i.3i33aa4a a *=a. 1694611 

c = o.36a5aoa5 , ; 

b*=o. 9371186 

a=i .67584367 
b=o. 95080317 


1 . 34*3335 
1.3433687 

+ 638 


r / *=o. 3913773 =[9.46445553] 
r**=o. 5983468 [9.73333776] 

0.889754® 

aV=o . 7583o3o6 

c*=o. i3i43og4=[9- 1 186646] 
c= [9.5593333] 


de là résulte la solution suivante ; 


=* [9 • 6762680] , r* = [9 . 88848 1 1 ] = o • 77353695 , 

d = [9-83738895]^ r 1 = [9.7322409] =0.53981000, 
c — [ 9 . 55 g 3 o 56 ] =0.3624980. 


5 
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a3. Voici, d’après celte solution, le calcul des ëlémens par 
les formules (14), (*5), etc. 

s — 0.8379336 3= [9.9233039], i— i"=o.oG/{ 3855 =fR. 80878803}, 
67708184]. J — / =0. 3981135=^9. 47458 oi5]. 


s — 9 47438 oi 5 

l/s 9.9616019 5 

A 9.4169831 

9.9118890 
A — B.. 9.058871 1 
cl / r ". . 9 . 60 354*1 1 5 

sin i 4-° 9.8553a49 5 
cos , 4 '°' 9 ■ 8434835 
cos’.... 9.6869670 

1° 9 . 88848 1 1 

n g. 0754481 


s-r°... 8.8087880 5 
V/11... 9 83 c 54 * ^7 7 

B 8.S473188 a 

9 . 43.79 8a 1 
9.3113467 a 

^ = o. 1636847 
o . 8373 1 53 


c 9. 55 g 3 o 56 

v//° g . 944 a 4°5 5 

c i / r ° 9 . 001546 1 5 

£(A-B)=I,A+/{i-|) 


Par les Formules directes, ftangj4'°=[o. 0118412], 
on trouve... I n =[9. 5754477] • 


Nous nous en tiendrons au premier" re'sul ta t, qui parait devoir 
être le plus exact; nous calculerons ensuite le temps employé 
à parcourir l'anomalie 4-‘, par la formule 


< = ( 3 no*(6H*+i® , H*) [..7644,79]; 


il en re'sulte t = 3 6.558q5 

e'poque de la i re observation octobre 22.68488 

passage au périhélie..... novemb. i8.o4353. 

Suit le calcul des latitudes et longitudes héliocentriques. 


<r* 

9.6762680 

a ' 

9.8575889 5 


0. 1 1 >5189 
9.6278110 


0.2677591 

9.3390147 

lang b ° 

lang b " 

sin£° 

9.4155979 

sia €' 

9.434*627 5 

cos £’ 

9.9847557 8 

COsê ' 

9.9833416 4 

long €‘ 

9.4508321 2 

tangf' 

9.4508211 1 

cot£'sin(a* — A’) 

0.2289950 

col i'sin (#' — A') 

0.3G47019 

sin(*° — A’). . . 

9.6598371 

si 11 ( a . ' — A').. . . 

9.8i55a5o 1 

a ° — A \ ... = 

37' 11' 18" 

Il 

* 

1 

40’ 5 o' i 5"45 

A° . — 

29. 19.50 
56 . 5 i. 8 

A' — 

41.22.26 

a° — 

a' — 

82. 12.41 .45 
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Pour calculer la longitude N du nœud ascendant, on a la formule 

‘ang^— N ) = ft) ta "g ~~T~ » dou rcsul,e 

— N = i 4 ° 52 ' 36 "; mais il est bon de vérifier celte valeur par 

la double équation col l = cot£°sin(a’ — N)=cotÊ'siu(a' — N); 
voici le calcul pour cet objet, en supposant 

— N= i 4 ° 5 a' 36 ''-f-jo"(.r), io''(x) étant une correction indé- 
terminée : 

*“ — N= 7 i*a 3 ' 44 ' , + I ° ,, (- r )) *' — X=97°5' i7"45-(-io''(x), 


sin(a° — N) 9. 9766909 -f~7ox, sin(œ' — N) 9.99666849 — a6x 
cotS* 0.56915788, cot£'.... 0.54917889 

coll o. 5458487 84-70X, cotl 0.54584758 — aGx; 

égalant ces deux valeurs, ou a i4-f-g6x=o, x= — 0.146, et 
io''(x)= — i’' 46; donc * 


— N= i 4 ° 5 a' 34' 1 54 j cot I = [o . 54584776] , 

N = 345. 7.25.46, 1= i5«53' o"54; 

enGn nous avons à calculer la formule cos II = cos cos (a* — N), 
pour avoir le lieu du périhélie : 

a’ — N = 7i°35'4a"54 

H = 72° 3 ' 3 2 "g 
cos(a*— N). . . . g. 5 o 38445 • 4 ” = 91.33.43.3 


cos£°... 9.9847658 N = 345 . 7.35.5 

cosll 9.4886003 lieu du pe'r. 148.44.4' -7 


En rassemblant tous ces résultats, on a les élémens de l’orbite 
comme il suit : 


Log n 

Passage au périhélie 

Inclinaison 

I.ieu du nœud asc.. 
Lieu du périhélie. . 


9.5754481 
nov. i 8 .o 4383 

i 5 * 53 ' o ''54 
545 . 7 - 25.5 
148.44.41.7 






( 36 ) \». 

a4- Ces éle'mens diffèrent asse* nolableinent des élémens. cor- 
rigés qu’on trouve pag. 3o du premier supplément; cependant 
on va voir qu’ils représentent plus exactement les trois observa-^ 
lions d’où ils sont tirés. 


Calcul du lieu du 3o octobre. 


Époq. donnée, oct. 30.68670 
Pass. au pér. nov. i8.o4 385 
t — 18.35713 

cos 14- 9.9829734 5 

cos* ........ 9 . 7859469 

n 9 . 57544 s» 

r 9,789501a 

sin I. 9 . 437246a 

•inK 9.9 991 57 4 1 

sinî. g. 4364o36 1 

co sC 9.9831620 5 

T 9.78 9501a 

SQ 9.7726652 5 

cos S. ..... . 9.9181410 a 

SI 9.6908042 7 

R. ........ . 9.9964960 


1 


V 

1 . a638o47 & 
n*. . 9.3631781 5 
T. . . 1 . goo63a6 3 
T = 79.548617 
4 = 77* 11' ao' 66 
pér.— Q = 163.37. 16. ao 
K = "86725755:54 
cos 1 9 . 9830989 * — N = 86.17.27.03 


a 

9.6g43o8a 7 TI 

0.4946617 

o.5o53383 

r. . . . * 

9.7895013 

sin C. 

g-4364o36 1 

cos T 

9. qa 101 12 a 

TI 

9.1469160 3 
9.700078a 3 

tang b 

9.4468 >78 

b, tsz 

1 5° 37' 5a* 6 

lat. observ. . . 

i 5.37. ai 

erreur... 

-+- 3i .6 

l’erreur est donc de 7", 6 suri 

Par les élémens corrigés, c 

33" sur la première et de 4 


cotK.... 8.7948540 
cot(» — N) 8.8117601 

voyez le triangle SQT 
dans la 6g. 3. 

R 9-9964960 

a 9.705 5822 5 


N = 345. 7 25. 5o 
• = 7 1.24. 5a. 53 
A = 37.19.48 
S = "34. 5. 4-53 

SQ.. 9.7726633 5 
sin S 9.7485107 9 
QI. . 9.52iit4cT 4 
TI. . 9 . 7000782 3 
9.8210958 1 


tangT 

T = 33° 3i' 8*4 
180° + A = 317-19. 48 
a = 183.48.39.6 
long, observ. . . 1 83 . 48 . 3a 
erreur. ... -j- 7 .é 
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a 5 . Quant aux lieux des aa octobre et 5 novembre, ils sont 
représenlés avec l’erreur de quelques secondes seulement par 
les élémens corrigés; mais, d’après la nature de notre méthode, 
ils doivent être représentés exactement par les élémens que nous 
venons de trouver. C’est ce que nous allons vérifier par des 
calculs qu’on doit regarder comme inutiles, mais qui pourront 
servir d’exemple à des calculateurs peu exercés. 

Calcul du lieu du a a octobre. 


Kpoque.... oct. sa. 68488 
Pass. aupér. nov. i 8 .o 4385 
t = a 6 . 358 g 5 

cos i 4 -. . . . . g. 8434834 3 

cos* g.G 86 g 668 6 

n q. 5764481 

r. g. 888481a 4 

cos S 9.9847557 8 

SQ 9.8733770 a 

cos S g.g 4 gi 5 c 6 6 

SI 9.8333876 6 

R. 9 . 9974 3 1 o 

9.8349666 6 
o. 6683935 5 
a 0.3317074 5 


t. . . 1 .4309381 
n*. . 9. 363 1731 5 
T... a. 0577559 5 

T = 114.233618 

»in I 9.4373463 
sin K 9.9783516 3 
sin C 9.4150978 3 


R-- 9 - 99743 '° 
a . . . 9.530755a 6 
TI... 9.818176a 6 


rsinf g. 3040790 7 

cosT 9.841906 6 1 

9.1459876 8 

TI 9.518 1763 6 

la n g b 9.6378114 a 

b = 33*59’ 53 *i 

observ 33.5g. 53 

erreur.. -J- 0.1 


4 = 91 ° 33 ’ 43 * 4 

r = 163.37.16. a 

K. ... = 7a. 3.3a! 8 


Cfisl. . .... 

tang K... . 

9 . 9830939 
0 . 4897507 

tang(«— N) 

0.4738446 

«—N. . = 

7 I*a3’ 4a” 43 

N — 

345. 7.a5.5o 

« ~ 

56 . 3 1 . 7 . 9^ 

A = 

«9 . i q . 5o 

S — 

37.H.17.93 

s<? 

9.8732370 a 

sin S * 

9 . 6698368 1 

Qi 

9 . 5330738 3 

TI 

9.518176a 6 

tang T. . . . 

0.0148975 7 

T — 

45*58’ 57' 

i8o°+A =s 

aoq. 19. 5o 

a . . . . . = 

i63.ao.53 

obsen-. . . . 

i63. ao.53 

erreur. , . . 

O 
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Calcul du lieu du 3 novembre. 


Époque. .. Nov. 3.70088 

Périhélie i 8 .c 4383 

t = 14.33390 


ï+ = 

33.24* 2. 6 

COsi 4 ' 

g.gai6o3G g 

COS* 

g. 8433073 8 

n 

g. 5754481 

r » . . 

9.7322407 a 

cos î. 

9. q8334 16 6 

SQ 

9.7155833 8 

cos S 

9.8788468 8 

SI 

9-5944=93 •> 

n 

q.gg6o53o 


9.5985763 6 


o.3<)66aiS 

a 

0.6033785 

rsinC 

g . 1 664o3a 6 

coa T 

g.g3ga54o a 


g. 105657a 8 

TI 

g. 7766428 0 

taog b 

g. 3390144 6 

b — 

13» 3' ag* 

ohserv. .... 

ia. a ag 

diff.... 

0 


t. . . 1 .1560334 g 

* - 

66*48' 5 * 3 

n*.. g. 363 i 7 ai 5 

r 

163.37.16. 3 

T... 1.7908603 4 

K.. ... = 

96.4g. 11 

T = 62.066934 



sinl.. 9.4373463 

cos I 

9.9830939 

sinK. 9.99691 63 4 

— tangK.. 

9.9333983 4 

sia C 9.4341606 4 

-taiig(a-IV) 

0.9053931 4 


N. . = 

97 ° 5 ' * 5.84 


N — 

345. 7.s5.5o 


« = 

8a. ta. 41 .34 


A — 

4 l .33 36 

R... 9.9960030 

S — 

40. 5 o. i 5.34 

a. . .. 9.7805898 a 

SQ 

9.7155823 8 

TI... 9.7766428 2 

sia S 

9.8155337 3 


Qi 

9-53i io 5 i 1 


TI 

9.7766438 3 


tang T. . . . 

9.7544633 9 


T — 

39° 36 ' io" 9 


i8o* + A = 

331 .33.36 


a = 

191.46.15.1 


observ .... 

igi - 46 .i 5 


erreur. . + 0.1 


a6. Il résulte de ces calculs qu’en effet l’orbite trouvée par 
la seconde méthode, satisfait aux trois observations données 
plus exactement que l’orbite corrigée de la pag. 3o, 1" suppl. 
Nous n'en conclurons pas qu'elle doit être préférée à celle-ci, 
car la véritable orbite est celle qui satisfait le mieux à l'ensemble 
de toutes les observations, et non pas seulement à trois obser- 
vations déterminées. 

Dans les cas ordinaires, il n’est pas nécessaire de pousser la 
précision aussi loin que nous l'avons fait dans cet exemple; 
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alors les calculs deviennent plus courts et plus faciles. A plus 
forte raison n*a-t-on pas besoin, dans la pratique, de rechercher 
une parabole plus approchée que celle à laquelle nous sommes 
parvenus. Cependant, comme cette question peut être de quel- 
qu’intérèl pour les analystes, nous allons faire voir quels moyens 
il faudrait employer pour la résoudre. 

On a du voir, par les calculs précédens, qu'il suffit de con- 
naître les deux distances accourcies a° et a', pour être en état 
de déterminer les élémens de l’orbite parabolique par un calcul 
rigoureux et purement trigonomélrique; car, par ces deux 
distances et les données de l'observation, jointes à celles qne 
fournit le mouvement de la Terre, on connaît immédiatement 
les deux rayons vecteurs r% r', et la corde c qui joint leurs 
extrémités. Ensuite des quantités r”, c, on déduit facilement 
la distance périhélie n , l’anomalie >{.* qui répond au premier 
rayon vecteur r°, et l’instant du passage au périhélie. Enfin les 
mêmes quantités jointes aux deux premières a’, a', font con- 
naître le lieu du nœud N , l’inclinaison I , et l’angle H qui donne 
le lieu du périhélie sur l’orbite H -f- *vj.' ■+■ N. 

Supposons maintenant qu’à la place des valeurs connues j% 
ff', nous mdtïom et x et y étant des cor- 

rections très petites dont on négligera les puissances et les 
produits au-delà du premier degré ; il en résultera des valeurs 
de r°, r, c, affectées des corrections r el/j ces valeurs étant 
substituées dans l’cqualion des temps, on aura une première 
équation entre x et y , au moyen de laquelle il ne restera plus 
à déterminer que la correction x-. 

D’après les valeurs de /•% r, affectées chacune d’une cor- 
rection de la forme kx, on déterminera les élémens de l’orbite; 
et par ces élémens, on calculera ensuite le lieu de la comète à 
l’instant de la seconde observation, ôn trouvera ainsi une er- 
reur sur la longitude qui sera de la forme 7 ", 6 -f -rnx, et une 
sur la latitude qui sera 5 i ',6 + nx;[il sera donc facile de déter- 
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miner la valeur de x , telle que ces erreurs soient anc'anlies en 
tout ou en partie, ce qui fera connaître la parabole la plus ap- 
prochée de la véritable orbite. 

27. Nous allons appliquer ce procédé à notre exemple. Soit 
pour cet effet dT°=m'x , cfo ’ = a'jr. 

' L’équation entre <r et r étant de la forme r* = A — Bs’ + Ca'*, 
on en tire 2rdr=(aCr — B) dr, ce qui donne 

dr° = — x(o.o8oai) , «/(/■**) = — x(o. 12/408), 

dr' = — 7(0 .18170), </(/*) = — 7(0 - 19617). 

Dans cette même hypothèse où l’on fait dr°=a°x et dr'=ir'j , 
on aura 

dV =x(o . 06753) — 7(0 . 335 a a) , 

^(c‘)=^* 4 ./‘- aV)=— ^(o.a 59 i 4 )-fX°- 474 2 7 )» 

* dc= — x(o . 35746) -f- 7(0 .65417); 

or, d’après l’équation des temps, il faut que la quantité. ... .-1 

1 1 

. (/•“-f-r'-f-c)* — — c)», soit indépendante des correc- 

tions x et 7 ; ainsi on doit avoir 

s i {dr* -f- dd -\-dc) — u 2 (dr° + dr' — de ) = o , 

ou 

£ 

o=—x(a . 43766)4-7(0 . 47 2 47 )— (“)’ Wo. a 77 a4)-7(o. 8358 7 )] ; 
substituant la valeur = [9.8769688], on trouve 
7 = x[g. 768343]. 

Au moyen de ce rapport, toutes les corrections s'expriment 
par la seule indéterminée x , et on trouve 

d = 0.7735370 — x(o.o8oai), 
d — o . 5398*1 00 — x(o . 1 o 6585 ) , 
c = o.36a498o 4“ x(o.oa 6 ag 5 ). 

Puisque nous n’avons plus qu’une seule indéterminée x, et que 
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le calcul des e'ie'mcns a déjà été fail dans la supposition de x=o, 
on pourra donner une valeur déterminée à x, ce qui simplifiera 
les calculs, et cependant les résultats étant comparés à ceux 
que donne l'hypothèse x=o , on en pourra déduire ceux qui 
répondent à une valeur quelconque de x. 

28. Soit doncx= loooounite's décimales du septième ordre, 
on aura 


r°= 0.7734568, d’où résultent les log n 

r' =3 o. 53g7o34 élémens de l'orbite périhélie... 

c = o. 36a5a43 qu’on joint ici. inclinaison. 

6767023] lieuduQ.. 

r =[g . 83782023 lieu du pér. 


9 . 57536809 
no v. i8.o3g6ga 

îS-SS'Sa"^ 
345. 7. 7.25 
148.47. 7.95 


Calculant, d’après ces élémens, la longitude et la latitude ponr 
l'instant de la seconde observation, on trouve l’erreur de la 
longitude — o",7 , et celle de la latitude 34 "; donc pour une 
valeur quelconque de x 

l’erreur sur la longitude = 7 "6 — 8",3 Ç~r~)> 

et l’erreur sur la latitude = 3 i", 6 -f-a ",4 (“^)- 

On voit qu'il n’est pas possible de faire disparaître les deux 
erreurs à la fois; d’ailleurs il convient de ne pas faire x trop 
grand, pour rester dans l'hypothèse des corrections très petites. 
C’est pourquoi on peut se bdrner à rendre nulle l'erreur sur 
la longitude, ce qui donnera x = 9igo, et alors l’erreur sur 
la latitude sera 35 ", 8. Ou trouverait ensuite facilement les élé- 
mens qui ont lieu dans ce système; mais ce résultat n’améliore 
pas sensiblement la solution du problème. 

28. En réfléchissant sur la seconde méthode que nous venons 
d’exposer et sur l'application que nous en avons faite à un exemple 
déjà connu , on voit que cette méthode présente trois degrés 
successifs d’approximation, qui s’obtiennent en donnant aux 
quantités A, A', e , des valeurs de plus eu plus approchées. 

6 
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Pour avoir la première approximation, on néglige les quan- 

8 V 

iilés du second ordre, et on prend simplement X' = g-^> 

e = o; alors soit qu’on emploie les équations (5) ou les équations 
(6), on trouve aisément, après quelques essais, une première so- 
lution qui sera le plus souvent fort imparfaite, mais qui servira à 
diriger les hypothèses dans le calcul delà seconde approximation. 

Pour obtenir celle-ci, on a égard aux quantités du second 
ordre; on prend donc les valeurs de A, X', e, comme il suit: 



on laisse alors r indéterminé dans les valeurs de «•’ et a', les- 
quelles se trouvent exprimées en fonctions de a et de r; chaque 
hypothèse faite sur la valeur de er, donne immédiatement celle 
de r ; on en déduit ainsi des valeurs entièrement connues de (r* 
et ff', lesquelles servent à calculer les trois quantités r°, r 1 , c, 
et celles-ci étant substituées dans l’équation des temps , on a l’er- 
reur de l’hypothèse. Enfin par deux erreurs semblables, suppo- 
sées l’une et l’autre assez petites, on trouve la vraie solution. 

Tel est ce que nous appelons le second degré d’approximation; 
il répond à celui que donne notre première méthode, simpli- 
fiée dans le premier supplément , et en général ce degré est 
très suffisant pour les besoins de la pratique. 

ag. La méthode de M. Ivory., qui coïncide avec celle de 
M. Olbers dans les cas ordinaires, donnerait un pareil degré 
d’approximation , et s’accorderait presqu'entièrement avec la 
nôtre , si l’auteur n’avait pas négligé les termes du second ordre 
dans la valeur des quantités analogues à celles que nous appelons 
X et X'. Ainsi M. Ivory suppose, comme M. Olbers, que la 
corde c est toujours coupée exactement en raison des temps 
8 et 8', par le rayon vecteur intermédiaire SC; or, cette sup- 
position , qui est permise lorsque les temps 8 et 8' sont peu 
différons l’un de l’autre, ne l’est plus, lorsqu’il y a une inégalité 
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sensible entre eux, et surtout lorsqu’en même temps les rayons 
vecteurs r*, r, sont plus petits que la distance moyenne de 
la Terre au Soleil, comme cela arrive dans l’exemple que nous 
avons calcule'. 

La méthode de M. Ivory dont nous venons de parler, est 
fonde'c sur le rapport des distances accourcies <r°, 7' , qu’on 
déduirait des équations (5) , en mettant R et A au lieu de 
Ri et Ai qui sont plus approchés ; elle a d’ailleurs l’inconvénient 
d’employer des auxiliaires inutiles, comme l’a déjà remarqué 
M. Burckhardt, dans la Connaissance des Tems de 181g. 

M. Ivory a modifié sa méthode pour la rendre applicable aux 
cas d’exception qui ont lieu , lorsque le lieu du Soleil est à 
peu près dans le même grand cercle que les lieux apparens de 
la comète dans la première et la troisième observation. Alors 
cette méthode se rapproche de celle que nous suivons, en fai- 
sant usage des équations (6); mais elle est encore compliquée 
par des auxiliaires inutiles, et elle suppose toujours que la corda 
est partagée exactement en raison des temps, ce qui doit l’em- 
pêcher, dans bien des cas, d’atteindre la même précision quo 
donne notre second degré d’approximation. 

3o. Enfin, dans ce que nous appelons le troisième degré 
d’approximation, nous avons pour but d’obtenir une parabole 
qui satisfasse .aussi exactement qu’il est possible aux trois ob- 
servations données. 

Pour cela, il faut calculer de nouveau les quantités A, A', e , 
au moyen des valeurs de r°, r, r', trouvées par le résultat de la 
seconde approximation , et à l’aide de deux hypothèses sur la 
valeur de 7 , on obtient aisément la solution cherchée. 

Les élémens déduits de cette solution satisferont toujours 
exactement à la première et à la troisième observation ; il restera 
donc à calculer par leur moyeu, le lieu gc'ocentriquc de la co- 
mète, à l’instant de la seconde observation, et ce lieu différera 
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toujours très peu du lieu observé, à moins que la véritable or- 
bite ne s'écarte sensiblement de la figure parabolique. 

Enfin nous avons fait voir comment on peut introduire, dans 
les élémens de l'orbite, une correction variable qui permettra de 
satisfaire encore plus exactement, s’il est possible, à la seconde 
observation. Cette recherche est au reste de pure curiosité, car 
lorsqu’on a obtenu des élémens qui satisfont exactement à deux 
des trois observations, et qui représentent à moins d’une minute 
la longitude et la latitude de la troisième, il est inutile de pous- 
ser plus loin. l’approximation. 

Application à la seconde comète de i8o5. 


3i. Dans le tableau de l’art. 35, 1 er supplément, nous pren- 
drons les trois observations suivantes: 


Époques. 

Longitudes. 

Lat. bor. 

Long. © — 18c”. 

Lng R. 

JÎOY. « 8.45203 
93.3294l 
Déc. 5.2g58i 

90® a 3 ' 37’ 
94.4,. 4 
a. 7.11 

3o° a' 40' 
97.95. 35 
3.90.45 

54” 1 1' 36' 
Si. 8.3o 
73.17. 4 

9.9947847 

9.9949071 

9.9933649 


Si au lieu de l’observation du 5 décembre, nous eussions 
choisi celle du 5o novembre , l’intervalle de temps de la pre- 
mière à la troisième eût été moins grand et moins inégalement 
partagé par la seconde observation , ce qui aurait produit sans 
doute des résultats plus exacts. Mais voulant mettre la seconde 
méthode à l’épreuve, nous avons préféré à dessein les circons- 
tances qui lui sont le plus défavorables , et qui au reste le se- 
raient également aux autres méthodes avec lesquelles on voudrait 
lia comparer. 

Nous allons encore donner le détail des diverses parties de 
l’opération , en faisant les calculs avec un degré de précision, 
que n’exigent pas les applications ordinaires.. 
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Dans cel exemple, on a 6 = 6.87038, 6' =11.97340, 
8 + 8'= 18.84378; ces valeurs exprimées en jours moyens doi- 
vent être multipliées par le facteur [8.a3558at], ce qui donnera 
les valeurs logarithmiques suivantes : 


6 9.0735638 

0 ' 9 - 5,5 799 6 


0+8'... 9.5107501 


r£r***‘ 9 - 56 i 8 ia 7 , 

rfp-"* 9- 8 °3°495. 


A ces premières valeurs, nous joindrons celles de Ai et Rr, 

données par les formules A 1 = j(A* + A + A') , 

logRi =^log (R°RR') , 


Ai = 6 a* 52 ' a 3 % sin Ai.... 9.949^892, 
log Ri = 9.9941 187 , cos Ai.... 9.6589501. 

Pour le calcul des quantités P et D, nous rassemblons ici plu- 
sieurs logarithmes constans qui auront leur usage dans les opé- 
rations ultérieures. 


tang 4 “ 9. 762:1170 
tang 4 g. 7 i 5 n 3 i 
tang 4 ' 8.7668757 


— sin(a' — a) 9.5840221 
»in(a° — a') 9.6610973 
— sin(o — a e ) 8 . g 143324 


sîn(Ai— a') 9.9407776 
»in(Ai — a ) g. 791165S 
sio(Ai — o*) 9.7416540 


D'après ces logarithmes, on peut voir, sans faire le calcul, 
que P est négatif; ensuite si on se borne à quatre décimales, 
011 aura en valeurs logarithmique» D= — [9 .346a] +[9. 376a) 
— [7.6813] = — 0.3319 + 0.2378 — 0.0048 = 0.0111; d’où 
l’on voit que DP est certainement négatif, et qu’ainsion a r>R; 
mais de là rien n'est à conclure pour la valeur de a. 

3a. Dans cet exemple, il convient de faire usage des équa- 
tions (5), dont le dénominateur est suffisamment grand, et qui 
sont alors les plus simples. On en tire 


a‘ = [9 . 94 54456 ], a' = — A — [9.17891 80] ^ 
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cela posé, sans bous arrêter au premier degré d'approsimation , 
qui ne sert qu'à mettre sur la voie des premières hypothèses, 
nous passerons immédiatement au second, et pour cela nous 
appliquerons les formules (9) qui donnent 

log (^) = ». 1969505 — p [7 .G5G5707] , 
log (~) = 0.4381875 — p [7.8191869J ; 

nous aurons donc pour calculer fl* et J par le moyen de fl et 
de r, ces formules 

• j 

<s° = 9 [0.1423961] — p [7 .603016] , 
a 1 = fl [916171053] — p [6.998105]. 

Ayant fait d’abord une hypothèse sur la valeur de fl, on calculera 
r par l’équation 

r* = 0.9736755 4- <r [o. aoo 6555 ] + fl* [o.to 556 ai]; 

connaissant fl et r, on aura par les formules précédentes les va- 
leurs logarithmique* de <f et fl', ensuite on calculera r 5 , r', V 
et c, par les équations 

^* = 0.9762689 + fl* [0.253795»] 4 - fl°* [o. ia 53 a 8 o], 

= 0.96990534- fl' [g. 8 o 33 g 5 o] 4 - fl'* [0.0014618], 

V sa o . 9 1 g 56 1 o -f~ fl* [9 ■ 85 iog 5 g] H- e-fl'fo . 9650522] 

4- a' [9.78341 16], 
e’ ss r“ 4 - r'* — 2 V ; 

enfin pour vérifier l’hypothèse ou en connaître l’erreur, on a 
l’équation des temps 

C(G + G') = (r* 4 - r 4- c)ï — (r* 4- r> — c)' = 1 - 9449 l8a - 
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33 . Soit pour première hypothèse a 
oq aura 


o. 157 *=[9.1958997], 


r=[o. 0^91 776] , r'=[<). 335448 *]] , .=[8.813*961], 


r°*= 1 . 4*7 1 8634 = [o . 1 5448064] r /, =i .01 54 o 58 =[o . 00663964] 

1^=1.19464896 0.077*403* r°*=i .4*718634 o.co 33 ig 8 a 


r'=i .00767346 

»*+»•'= a . ao* 3 a * 44 . a^= 4 . *890*37 

0=0.43 7485*6 b î =a> 3445353 


a=a. 63980768 
b=i .76480716 


1 .9444874 

1 .9449183 


*.44*59*14 

aV=a.* 5 i 19878 

c’=o. igi4*i46=[9 *819*68] 
9.6409634 


erreur — 43 o 8 


Soit pour seconde hypothèse <r=o. i 5 ya =*[g. ig 645 a 5 ] , on 
aura 


r=[o . o 4 g* 4635 ] , ^=[9 . 336 ooa 3 ] , ' * =[8 . 81 * 84 g 3 ] , 

r“*=l .4*78460 =[o. l 54880 l 6 ] »^=l . 0 l 54 Sg* =[0- 00666676] rfj 

r*=i .1 g 4 ga 336 0.07734008 1.4*784*0 o.oo 333338 

JK 


r'=i .0077049* 


Lfc'i: 




a. 4433 i 


*0*6*898 • ï =^.9 S miS3;- ■ 

=0.4376681* b a . =a 34478 , 4 r*=o. 19. 55 i 4 =[g. *8**8988] 

,=*.640*9640 —7 r=r â *4>'4494 _ 


a =* . 640*9640 
1)=1 .76496016' 


1.9454® 9 

1.944918a 


erreur -+■ 5i57 


{> 

Otj I uVj.iriOl,, . ; 'T' Ü0I8 

vjèiîio'i 


P.Sj jc]oo.6 ] r ecîo* 1 r\ ‘ i-«r-'*V - 

delà résulte la solution suivante 

* ^*=[9. 3357004], r*= [0.0773856] , e= [9.64109935], 

<r =[o. iq 6 i 5 i 3 ], r =[0. 0493080] , 'HH 

<r'=[8. 8135476], r'=[o. oo 33 a 6 o]. 

"’ T l ' f '. ' . . {;. . I ‘ ^ V ' * 

Cette solution serait suffisante pour les cas ordinaires, parce 
que nous avons tenu compte des terme* du second ordre ; mais 
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nous pouvons obtenir une plus grande approximation, en tenant 
compte des termes du troisième ordre. 

54. Pour cet effet, il faut déduire du dernier résultat les va- 
leurs corrigées de A, A', e, d’après les formules 


£ = 


mV‘ 


6 (r/j 

J68' 

rrP» 


l 


mô* 

6 (rV)* 


log A = log , 

log a'= log ( 6 -^f) 6 qp s - » 


on trouvera aiusi les valeurs logarithmiques 


1 — e =[9. 996062a], A=[g. 565 oa 5 o] , A'=[9.8oa353g] , 

- . ' • * 

au moyeu desquelles ou aura 

0° = «■ [0.139x539], e 7 = «■ [9.61 ig 55 a], 


ou simplement, parce qu’on n’a plus besoin de <r, 


O* ¥= O* [ 9 . 47 a 8 oi 3 ]. 


( Nous négligeons ici la correction de Ai qui pourrait se dé- 
duire des équations de l’art. 10 ). 

Soit pour première hypothèse ff* = [9. 334 a 1 5 a]; on aura 
e' = [8. 80701 65 J. 


r"=i .415731 13 =([0.1540376] 
r*=i . ig4o3g83 0.0770188 


r , =i .00737557 


r*-f-r'= a.aoi4i54o a*=4- 18734*4 
0= 0.43766154 b i_ a . 3435775 


«=9.63907694 
b=i .76375386 


1 .9448649 
1.944918a 


etr. 


— 533 


r ■= i . 0 1 48 o 55 a = £ o . oo 638 a 8 a] 
1.42573113 0.00319141 

a.44o53665 
aY=a. 34898896 

c*=o . 1 g 1 5476g=[9 . 38237684] 

9. 6411 384 » 


V* 
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Soit, pour seconde hypothèse, a” = [g, 3343ooo] , on «tira 
o' = [8.8071013]. 


r°*=i ,4a583log =r[o . 1 54o68o6] 
r°=i . 1 9408 168 o . 07703403 


r'—i .0073803a 


r’+r'^a . ao 1 46300 a >^4 . 3874345 
0= 0^3768963 b J =a-34 , 4l43 


a= 3 . 63 q 1 5 1 6a 
b=i .76377338 


1 .943010a 
1 .944.918a 


» / *=i . o 1481609= [0.00 63 86gaJ 
1 .43583109 0.00319346 

a . 440646 1 8 
aV=a . 349074°° 

c ! =o. igi5 7 ai8=[g.a8233a44] 
9.64116833 


err. + 930 


de là résulte cette solution: 


a’ = [9.3342463], r» = [0.0770244], 
ff' = [8 . 8070476] , r' = [o . oo3 19217], 
c =[9.64114862]. 


35. Pour en déduire les élémens de l’orbite, on formera 
d’abord les nombres s, u, s—r •*, s — r', et leurs logarithmes, 
par lesquels on calculera 4° et n , comme il suit: 


r*= i.iq 4 o 55 aa 
/ = 1 .00737733 
c = 0.43767183 

as— 3.63910438 
s = 1.31955319 

u= o.88i88o36 

r 

•iug4° = — — - 


>. . . . 0.1304265 7 
u. ... 9.9454096 8 
su.... o. o65836a 5 
s — r°=o. 1354969 7 
s—r’ =o.3iai 7 48 6 

’s — (,t — r*) \/u 

cyi*. 


C 4 *)’ = 4-38 7 3og4 

(»“)* = a. 3433910 

t. 944g 184 

(s— t°) 9.098633a a 

0 “ O 9 - 4943: 9 79 G 

(s — r*)(j — /).. 8. 593o3 ïï "8 

« 9.6411486 a 

Vr * o.o385iaa 


cÿr° g . 6796608 a 


r 
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(♦-/•O 9-4945979 6 (*—»')•••• 9 - 09 86553 a 

^/f o. 060a 1 3 a 8 v/u 9,9737048 4 



9,55461 1a 4 

g. 07*3380 6 


9.8a6g5a3 4 

9.554611a 4 


g. 38 1 5635 8 

9.5167368 a 

c^/r" 

9.679G608 a 

o. 5 a 86448 6 

mnx 4 * 

9.7019037 6 

0.6713551 4 

cos f4°- . • 

9.9365441 7 

Par les formules directes , 

tangi 4 °--*- 

g. 7 G 53585 9 

on trouve 

COS*7 4 ° 

9.8700883 4 i an gi 4 * = [9.7653582 5 ] 

7 * 

9.0770344 n 


n 

9.9501137 4 


Nous prendrons par an milieu 


n = 

= [ 9 - 9 5o,ia 9 ]» 

^ 4 * = 3o“i3'a7"8, 

, Ungi4' = 

= [9.7663584], 

4 * = 6 o.a 6 . 55 .^| 

Pour avoir l'instant du passage 

au périhélie, nous calculerons 

directement la 

formule 


t = (an’)ï(tang t 4 4 " 5 l an g 5 ï 4 ) t 1 • 7 6 44 1 79 5 ] > 

(an J )î. . . . 

.. 0.0756843 5 

1)... = 40. 3 i 4444 


1.7644*79 5 

a). . . = 4 • 560960 

<an g. 4 * i • 

. . 9.7655584 

44-875404 

0 

.. 1.6054607 

nov. . . . i6.453o3 

tang* 

, . . g. 5307 168 ' 

dcc. ... 3 i .337434 

1 

J* • * * * * 

. . 9.5338787 5 

passage au périhélie. 

3 ) 

.. o. 65 go 56 a 5 
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Suit le calcul des latitudes et longitudes héliocentriques: 


•" 9.3343463 

i:r ° 9-9339756 

tang 6°.... 9.7633170 

9.0194389 

ce» C* 9-9976i30 5 

tang C °. ... 9.0318368 5 

cot 6 * 0.3377830 

sin(A° — a”) 9.6336778 5 

*in(A°— •*) 8.8833877 

A*— <*= 4 * 33' 34 ' 71 
A«= 54.11.36 

-*= 49 - 49 -Ü-S 9 


r 8 . 8070476 C°=6* o' 1 O* 305 

tir ' 9.99680783 C=o. 13.47.640 

rang b ’.... 8-7668707 C’+C'=6TTr5r‘845 

«inC 7.5707311 3 0* — C=5 .47.33. 565 

ce» 1, 9-.q999969 9 

tangC’.... 7.5707341 4 

cot 6' 1 . a 33 1 fl 43 

ain(A'— a ') 9.9760979 6 

»in(A' — «') 8.7799564 
A' — 3*3/14*97 
A — 7 3 | 7 - 4 
*•'= 69:49.49.03 
“’= 49 -4.9 * 1 • a .9 


ï («'+«*)= 59.49.30.16 
i(« — «•*)= 10. 0.18.87 


II faut maintenant calculer le lieu du nœud N par la formule 


tang[N — i(a'-f.ot*)] = 


»in (C* + f) 

sin (C° — T) 


tang 



tang ï C*' — **)•••• 9-a4655io 9 as= 10*43' 


siu (£*-+-€') 9.03454084 ï( a 4 " a ’) 5 g. 49 . 3 o.i 6 

8.3810919 3 N= 70. 3a. 5 a. 89 

sin ( 5 *— €’) 9.00578668, ou N=a5o.3a.5a.89 

tango: 9. 3773063 5 


La seconde valeur de N est celle qu’il convient d’adopter, parce 
que siu (œ* — N) doit être de même signe que coté*. Mais, comme 
cette valeur est déduite d’angles assez petits, nous allons la vérifier 
parla double équation cotl=cotff*sin(a° — N)=col£'sin(*' — N). 
Soit, pour cet effet, N = a 5 o* 3 a' 53 *+io"(x), io"(or) étant 


( 5a )• 

une correction indéterminée, on aura 

■° — N= i 5 g* îff 18*29— io'(x). N= 179 0 16* 56 *o 3 - xo'(x), 
ao.43.41 -71 + 10*(x), 0.43. 3.97+10 (x). 

sin(«»— N) 9.5489248 1 sin (« — N) 8.0978509 6 + 16810X 

*ot i° 9781731 5 cot C.... 2.4292^58 6 

cotl o. 527097g. 6 + 556 j cotl 0.5271168 2 + 16810X 


Égalant ccs deux valeurs, 011 a i 6 a 54 as=>— 188.6, x=— 0.01 16, 
io''(a , )= — o*« 16, et par conséquent 

N== a 5 o° 3 a' 5 a" 884 , 1 = i6 «32'48 ''oG5 , sinl=[g. 45453468 ] 


cotl— [0.5270973 a] 

Enfin nous avons à calculer 
cos H = cos G° cos («” — N) 

a 0 — N = >59" 16' i 8"4 

— cos (<t° — N) 9.9709369 3 

cos 9-9976120 5 

— cos H 9.9685489 7 


cosI=[9- 981 63 ao ] 
H par la formule 

H = i 58 ° 37' 24" 47 
•4* = 60.26.55.60 
N = 25 o. 5 a. 5 a .88 

lieu du pér. 109.37.12.95 


D’après ces résultats , voici les élémens de l'orbite .- 


Cog n •/•/.••• 

Passage au périhélie.. 
Inclinaison. . 

Lieu du nœud 

Lieu du périhélie.... 


9.9501 129 
déc. 3 i .327434 
16” 3 à' 48 v 
a 5 o. 3 a. 52.88 
109.27. 12.95 


36 . Ces élémens diffèrent assez peu des élémens corrigés ,' 
pag. 48 ; ils doivent satisfaire exactement aux doux observations 
extrêmes, niais il faut voir jusqu’à quel point de précision il» 
représentent l’observation moyenne. 


« 
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Calcul du lieu du a 3 novembre. 


Époq. .nov.a 3 . 3 fl 24 i 
Périh... déc. 3 i .327434 
t = 38 . 003024 


= 26 “ 5 o' 25*225 

cost’Î' 9-95o4954 1 
cos*. . 9.9009908 2 
n. ... 9.9501129 


r — lien du pér. — Q 
K = i — r 
Ni lieu du noeud dcsc. 


r. . . .. 0.0491220 8 

sinl.. 9 . 454534 S 8 cos I g.g 8 i 63 ao 

sînK " ^*4°65838 4 tangK 9.42118695 


t 1.5798409 8 

-3 

n *... 0.0748306 5 

T 1.6546716 3 

2 = 45.15144 

4- = 53°4o'5o*45 
r = 38.54.20.07 

K = 1 4 . 46T3ÔT3S 

Ni— « =5 14*11' 18*41 
Ni 70.39.52.88 


ainC.. 8.8611185 9 tang(Ni— «) 9.4028189 5 
cosC. . 9.9988514 5 

T 0. 0491220 ,£ 

SQ. . . 0.0479735 3 


» = 56 . 21 . 34.47 

A = 61 . 8 . 3 o 

s = 4.46T55753 ’ 


cos S.. 9.9984855 6 
SI.,.. 0.0464588 9 
9 - 994 =<° 7 i 

o. 0529517 9 

1 . 1278508 8 

r 0.0491220 8 

sin C. . 8.861 ll 85 2 
cosT.. 9.9055868 


8.8158274 

TI. . . g. 1009107 9 
tangi 9.7149166 1 


R 9 -9943071 

9.1067036 9 

TI... r . 9.1009107 9 


. • ■>- 

b = 97 * 94 ' 56 ’ 84 

observ. 27 . 25.35 

erreur. — 38 . 16 


SQ. . . 0.0479735 3 
sin S . . . 8 . 92C9908 9 

Ql. . . 8.9689644 9 

TI. . . - 9.1009107' 9 

«angT... g. 868 o 536 3 
T = 36 °a 5 ' 3 7 * 76 
6k 8.3o 


obs. . 


crr. .. 


24.43.02.24 

H 4 s- 4 
+ * - 48-34 


Ainsi les éle'mens trouvés, qui satisfont exactement aux deux 
observations extrêmes, ne représentent l’observation moyenne 
qu’avec une erreur de 1' 48"sur la longitude, et de 38 " sur la lali- 
tude. On voit néanmoins par le petit tableau de la pag. 47, que 
ces trois observations sont encore moins bien représentées par 


( 54 )• 

les élémens corriges; ainsi malgré toutes les difficultés que pré- 
sentait cet exemple, difficultés qui seront rarement réunies à un 
pareil degré, la seconde méthode nous a conduits directement à 
une solution égale, ou même supérieure en exactitude, à celle 
que nous avons trouvée eu combinant ensemble un plus grand 
nombre d’observations. Il semble qu’on ne peut rien désirer de 
plus dans la solution analytique du problème des comètes. 

Nous avous établi la solution précédente sur les équations (6), 
nous aurions pu aussi l’établir, et peut-être avec quelque avan- 
tage, sur les équations (7) et (8). Pour donner uu exemple de 
l’usage de ces dernières équations, nous supposerons qu’on est 
déjà parvenu à ce que nous appelons le second degré d’approxi- 
mation , c’est-à-dire à la solution de l’art. 53, et d’après cette solu- 
tion, nous déduirons des équations (7) et (8) le troisième degré 
d’approximation. 

Nous avons déjà trouvé A= [g.’ 563oa5o], A'=[9.8oa553g J, 
on a en même temps e= [7,95551 19], E=[8. 1039764], 
E — e= [7.56aa3i6], 1 — e=[g. 996062a] ; substituant ces va- 
leurs et celles des quantités connues dans les deux premières des 
équations (7) et (8), on aura 

j 397365 ] — [6.4579480]; 
a' = a [ 9 . 5768860 ] -f- [7 . 6926 1 82 ]. 

Soit pour première hypothèse ^=[9. ig5oooo] ; on aura 
^"=C9 . 334 1 5g36] , r =[8.806641 4] . 


r**=i .4*566533 =[0.15401756] /‘—t .oi4763ao=[o. 00636471] 



r°=t. 19401231 0.07700878 

1 .4*566533 

o.oo3i8a35 

» 

^=1 .00735456 

3.44042853 



r'+r' =0.20136687 a , =4. 2871913 

aV=2.a4?85698 



=0.43768890 b*— 2.34**266 

c*=o. 19157155= 

: [9.a8a33ioa] 
9.641 i655i 


«=2.63905577 ,.044,647 

b=i .76367797 




err. •+■ 46 5 
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Soit de nouveau »=[g. »9497°o]i on aura 

^=[9.33412933], r =[8.8oG6l37l], 


r**=i .4 s 563 ooi =[o.i54co68a] j / *=i .oi476oo7=[o.oo636338] 

r*=i . ig3gg755 0.07700341 i.4»563ooi 0.00316169 


i'—i .co7353oa 

. 80 135 C 57 
C= 0 . 43767863 

a= 3.63902920 

b=i .76367194 


a 

1 ^= 4. 3871263 

b*=3. 3423146 

*•944.9 "9 
». 94491 8a 


3 . 44°3gco8 
2Y=s. 2 {883748 

c*=o. igi56a6o=[g. 3833107] 
9.641 1 5335 


err. — 63 


De là résulte cette solution 

ff”.,.. 9.3341329 2, 0.0770040 5 

e'.... 8.8066170 3 , r' 0.0031817 7 

c 9.64* *565 7 


. Voici, d’après ces données, le calcul des élémens de l'orbite: 


** = *■ *9 3 .999 a 4 

/ SX 1.0073533 o 
c ss 0.4376798 7 
S 5 = 9 . 639 o 3 a 3 I 

j 71 .3196161 5 

u = o. 88 i 8369 8 


s o.iflo4»47 1 

u 9.9455879 6 

su o.c658os6 7 

(i— r°) ...... 9 . 0987099 6 

($— /) 9.49438 13 6 

(s — r°) (%— /) 8. 5900836 ~2 


1 

(3i)“ = 4-a87l34l 
(au)* = 3 .5433160 
1 .9449181 

e 9.6411 565 7 

l /r*. . . . o.o385oao a 
«V**... 9 . 6796685 “^ 


(1—0.. 9. 4g 438i 3 6 
y*. . . . 0.06113073 5 
9.5546887 1 

.t) 9.82G9138 3 

9 38 idoi 5 3 
c(/r°. . . 9.6 796585 9 
»in i4°-- 9.7018429 4 


(. — r*) 9.098703a 6 
\/ u . .. 9.973 6939 8 

9. 0713962 "4" 

9.5545 887 1 

9.5168075 3 
0.3387059 a 
x) o. 67 12940 8 


^ar les formule. 


Ung;4-“= 

n= 


1/(1— r-)V(,_7j 

(J/r— Vu)* * 


on trouve directement 


co. 14*.. 9.936664 4 4 co.’l 4-'’ 9 873 1288 8 

Ungi-}, 0 9.7653785 r* 0.0770040 5 

n 9.9601339 3 

par un milieu n. . = 9.9501329 6 


tang £■!».. 9.7652785 
n 9.9001329 9 

44» sa 3o«i3' 11*38 
4* — 60.36. 32. 56 


( 56 ) 

« = (on 1 )* (tangi+H- itang’J 4 ) [1.7644179 5 ] 


(an 5 ) 1 0.0757144 4 

1.7644179 5 

langi^ 0 .... 9.7659785 

1) i. 6 o 54 ic 8 g 

lang* 9.53 o5570 

9.5898787 5 

a) 0.6588466 4 


1) 4 ° • 3 cg 8 a 3 

a) _ 4 .558759 

t s= 44 . 86858 a 

DOT l6.43303 

déc 3 i , 3 ac 6 ia. 


Passage au périhélie. 



r . • , . • 

. 8.8066170 3 

P — 6* d 5*558 

>:r* 9-9 aa 99Ô9 5 

i-y... 

. 9.9968183 3 

S =0. 19.46.8q7 

langé» 0 ... 9.7693170 

tang b'. 

. 8.7668767 

P+P.. = 6 . 13 . 53.455 

lin C 0 .... 9.0193408 7 

•in C. . 

. 7.5703109 6 

P— = 5.47.18.661 

cosC*. . . . 9.9976130 8 

cosC. . 

■ 8 - 999997 ° 0 


tangï 0 ... 9.0917397 g 

lang C 

. 7.5703139 6 

sin (0°-+-^) 9 .o 344366 6 

coti 0 .... 0.3377830 

cot b'.. 

. 1 . a 33 1 243 

ein (S*- — ”) 9.0007055 9 

8in(A 0 -o*) 9.6296778 5 sin(A'-a') 9.9760979 6 

o.o3o73io 7 

6in(À°-«°) b. 8821936 4 sin(A'-« 

) 8.779536a 3 

*angj(«'-« 0 ) 9.3466041 8 

A" — «* = 4 ° 33' 3 i* 9 g 6 A' — «' = 

= 3 ’ 97' 3*936 

tangx. ... 9.9773353 5 

A° =54. 1 1 .36 

A' = 

=73.17. 4 

X — io* 43 ' a 5*34 

«* = 49 - 49 -i 4 - 7°4 

«' =69. 5 o. 1.074 

59.49.37.89 


-1° 

4 . 9 - 49 - > 4 - 7°4 

70.33. 3 .a 3 

1 

% 

(•+-")= 

=59.49.37.889 

N =a 5 o. 33 . 3 .a 3 


!(.' — *°)= 10. o . a 3 . 1 85 

Pour vérifier la valeur de N, soit N — 25o° 55' 3’'+ 10 " far) 

m ° = 49 * 49 ' ' 4*704 « = 6g* 5 o' 1*074 

N = a 5 o. 33 . 3 -f»io*(x) N — a 5 o. 33 . 34 -io*(x) 

N = 1 5 g . 16. 1 1 .704 — 10" (x) » — N = 179.16.58.074— io"(x) 

ao.43.48.996+io"(jr) 0.43. 1 ;ga6-{-io*(x) 

1111(11*— -N) 9.5489614 3 + 556 x ain(«' — N) 8.0978073 5 +i 683 ox 

cotî* 0. 978967a 1 cotT 3 .4996 860 4 

cotl 0.5373986 4 + 556 x oot I o. 597 ig 33 _ 9 +i 6839 x 


Egalant ces deux valeurs, on a x=o. 02164, et io''(x)=o*2i6; 
donc 


N = a5o°33' 3*ai6 1= i6*3a'3o*89 »inl = [9.454419g j 

eotl s= [0.5973998 4] N= 159.16.11.488 co» 1= [9.9816497 4] 
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—cm («• — 9.97093135 9 H = i 58 ° 37 ’ 19 ' 36 

coa ff" 9 .9976130 8 4 '“ — Go.aG.aa. 5 S 

— cosH g. 9685444 7 N — « 5 o . 33 . 3 .aa 

. lien du périhélie 109.26 . 45 .<4 

Voici donc les nouveaux élémens de l’orbite que nous mettons 
en comparaison avec ceux qu’a fournis la première solution. 


D apres ces seconds élémens, nous allons calculer le lien 
du a3 novembre. 

Époque... Nov. s 3 . 3 aa 4 i t . 

Périhélie . . Déc. 3 i . 3 ao 6 1 a ■ 

* = ^7- 99820a" T, 

H aG‘ 5 o' 7*37 r=Pér. — g (N.) T, 

coeiih 9 .g 5 o 5 i 44 1 N 1= 70“ 33 ' 3 * aa 

co»*.. 9.90102883 K («rg. «le l*t.)= 4 ~-r r 

n 9.q5oi3aq 6 K 

r. 0.0491041 3 

»inl. 9.4544139 cosl 9.98164374 N 

•iuK 9.4066033 3 tangK.. 9.4310077 8 N 

•inf 8. 861016a 3 tang(Ni-«) 9 . 4 oa 85 o 5 a ». 

coeC. 9.9988519 9 A, 

r 0.049104» 3 R 9.9943071 S. 


1.5797630 S 

0.0748005 6 
1 .6545636 1 
45.i4oai4 
53 “ 4 °' > 4*74 
38 . 53 . 41 . 9 » 

14.46. 3 a. 8a 


9 ■ 998 4867 a TI. . . . 

0.0464438 4 r 

9-99 < 3 °7» * in f- • ■ 
o.o5aa357 4 cos T. 
1 . 1378094 8 
37“ aâ' 4* TI.... 

37.35 35 tang b 


observ 


. ' '.1 « ; 

Première solution. 

Seconde solution. 

Log n 

Passage aa périhélie. 

Inclinaison 

Lieu du nœud 

Lieu du périhélie . . . 

9.9501139 
déc. 5 i .337434 

16“ 3 a' 48' 
a 5 o. 3 a. 5 a .88 
109 .37. 13.95 

9.95013396 
déc. 3 i. 3 ao 6 ia 

iS“ 3 a' 3 o’ 8 g 
a 5 o. 33 . 3 ,aa 
109. a 6 . 45 . 14 
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Ainsi dans celle seconde solution l'erreur augmente de n'^Sur 
la longitude , et diminue de 7" sur la latitude. 

4o. Nous avons déjà fait voir que les élémens de l’ctrbile peu- 
vent toujours se déduire d'une manière rigoureuse, des deux 
distances accourcies s°, a', supposées connues j nous avons prouvé 
de pins qu’on peut ajouter une correction indéterminée aux va- 
leurs de ces deux distances, mais que pouè satisfaire aux équa- 
tions du problème, et particulièrement à l'équation des temps, 
il faudra que les deux corrections ne dépendent que d’une seule 
indéterminée. Cela posé , si d’un système déterminé a°—p f 
a' =7, qui satisfait à toutes les équations du problème, on dé- 
duit les élémens de l'orbite , et que par ces élcmens on calcule le 
lieu de la comctc à l’instant de l'observation moyenne , suppo- 
. sous qu’on trouve l’erreur E sur la longitude, et l’prreur c sur la 
latitude. Si d’un second système a , —p-\-<5'p i {-Stf, on 

déduit semblablement de nouveaux élémens, lesquels donnent 
pour la longitude l’erreur E+JE, et pour la latitude l'erreur 
c + eTc, on pourra en conclure que pour un système quelconque 
peu différent de ceux-là, dans lequel on aurait e 0 = p -f- x£p , 
<r'=</+xj'q, les équations du problème seront encore satis- 
faites , et que les élémens de l’orbite qui en seront déduits , don- 
neront l’erreur E + ar«TE sur la longitude et l'erreur e+Xtfe sur 
la latitude ; et celte conséquence a lieu pour toute valeur de x po- 
sitive ou négative , pourvu quelle n’excède pas une certaine limite 
où les valeurs simplement linéaires des corrections cesseraient 
d’être suffisamment exactes. . 

On voit par conséquent que par le moyen des deux orbites 
que nous avons déterminées , on' peut eu trouver une infinité 
d'autres peu différentes de celles-là, dans lesquelles les erreurs 
sur la longitude et la latitude du lieu moyen , seront représentées 
parE+jr<fE et -e-t-xje. Il s’agira ensuite de choisir x, de manière 
que ces deux erreurs soient atténuées autant qu’il est possible, 
pour obtenir la solution la plus exacte du problème proposé. 
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Nous avons trouvé dans la première solution E=io8",a4, 
t= — 58', 16 , et dans la seconde E+<fE=i i8’,g6, e+Sex= — 5i% 
de là résulte un troisième siytème dans lequel l’erreur sur le lieu 
moyen sera 

en longitude = io 8 ',a 4 + >r ( 10 *> 7 a )» 
et en latitude. = — 58",t6.-f-#( 7 ",iG). 


S’il y avait une valeur de x , peu éloignée des limites x—o, .t=i , 
qui rendît nulles ces deux erreurs, on aurait par les élémens 
eprrespondans, une parabole qui satisferait exactement aux trois 
observations données; mais on voit immédiatement que cette 
condition est impossible à remplir, car on ne peut diminuer une 
des erreurs sans augmenter l’autre. 

Suivant la méthode des moindres carrés on aurait à très peu 
près x= — îj, ce qui donne l'errcnr i'ag',5 sur la longitude, 
et — 5 o ',7 sur la latitude. C’est le résultat qu’il parait convenable 
d’adopter. 

4i. Quelle que soit au reste la valeur que l’on prend pour x , 
pourvu qu’elle soit peu éloignée des limites x=o , x=t , les 
élémens de l’orbite correspondante se trouveront très facilement 
par le moyeu du tableau précédeut, qui offre le résultat des deux 
solutions. En effet, si un clément a pour valeur A dans le pre- 
mier système, et A+JA dans le second, sa valeur correspon- 
dante à une valeur donnée de x sera A-f-arJA. Ainsi, en fai- 
sant x — — i j conformément à la détermination précédente, on 
aura les élémens suivans de l’orbite. 

/ 


I'°g n 

Passage au périhélie. 

Inclinaison 

J.ieu du nœud. ..... 

Lieu du périhélie. . . 


g.g5oo778 
déc. 5 i. 53 gÎ 7 
i6*33' 18 ' 
a5o. 3a . 55 
iog.a 8 . a 


Si l'on calcule d'après ces élémens les lieux de la comète pour 
les époques des trois observations données , on trouvera les diffé- 


(6o) 

Te n ces suivantes qui portent le signe -J- lorsque le résultat da- 
calcul surpasse celui de l'observation , et le signe — dans le cas 
contraire. 



i r * Observ. 

a* Observât 

3* Obatrvat. 

Longitude 

- «"4 

-t- i' 33" 6 

— o"8 

Latitude 

O 

— 5i-7 

+ 2.5 


Suivant la théorie précédente, les erreurs devraient être nulles, 
pour la première et la troisième observation, et relativement à la 
seconde, elles devraient être i'ag',5 et — 5o",7; si ces conditions 
ne sont pas exactement remplies, c’est que la valeur suppo- 
sée .r= — 1 | est assez éloignée des limites ,r=o,a:=i, pour que 
les termes du second ordre qu’on a négligés, produisent quelques- 
secondes d’erreur dans les résultats. 

Remarque. 

On peut supposer que les trois observations qui servent à 
déterminer l'orbite ont été corrigées de l’aberration et de la nuta- 
tion; mais pour obtenir la solution la plus exacte qu’il est pos- 
sible, il faut encore qu’elles soient corrigées de la parallaxe , ce 
qui ne peut se faire que lorsque les distances de la comète à la 
terre sont à peu près connues, c’est-à-dire, lorsqu’on est parvenu 
par la méthode précédente au second degré d’approximation , ou 
au moins au premier. Ainsi avant de procéder aux derniers eal- 
culs d’approximation, il faudra, au moyen des résultats précé- 
dens, corriger les observations données des effets de la parallaxe , 
et rectifier en conséquence les coefliciens des équations à ré- 
soudre , ce qui sera un assez grand surcroît de travail. 

Ou pourra cependant négliger cette correction, lorsque la 
distance de la comète à la terre sera assez grande pour que sa 
parallaxe horizontale ne s'élève pas à i', et dans tous les cas. 
aussi où l’orbite déterminée par notre méthode ne sera cousi- 
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dérée que comme une première approximation pour parvenir à 
l'orbite qui satisfait aussi exactement qu'il est possible à l’en- 
semble des observations. 


TROISIÈME MÉTHODE.. 

4a. Il nous reste maintenant à discuter une troisième méthode 
dans laquelle on suppose connues la longitude a et la latitude b 
de la comète, observées à un instant déterminé , avec leurs coef- 
ficiens différentiels du premier et du second degré, pris par 

. da db dda ddb - 

rapport au temps, savorr ^ , —j--. La question consiste 

donc à déterminer les éfémens de l'orbite réelle par le moyen 
d'une portion infiniment petite de l'orbite apparente, supposée 
parfaitement connue. Ce problème peut se résoudre dune ma- 
nière rigoureuse et fort simple par les formules de notre première 
méthode, dans lesquelles on supposera 9 et 6' infiniment petits; 
on pourra même sans inconvénient faire 9' = 9. 

Dans cette hypothèse on aura d’abord f =/■-}- F'8 — F"8* , 

S'=S-f-G'8— G'8*, et puisque /= lang a , g — , on ea 

déduira immédiatement i 


F' 


F" 


df i da 

d{ co b* a * dt 9 

i dda sina /da\* 
scos'a * dt * 2 C 03 s a\dt / 9 



sinatang& 

cos*a 


da . 1 db 

' <fi * cofàcosa * dt r 


1 ddb sinfl dadb sioA 
ncosacos'b * dt * cos*acos*b * dt* cosacos { b \</f / 

sinotangi dda tang&(i -f-sinVi) /da\* 

aco s’a " dt “ 2C0‘r , a \3îJ * 


D'ailleurs on a aussi F =f, G=g ; et comme il faut supposer 



C 6a ) 

dans nos formules < 1 = 0 , il en résulte 


_ r?/ da T?ti i ddû 

— 0 > ^ ~ ’ fa > * * • rit* > 


dt‘ 


db 




G": 


i Mb 
a cos* 6 ’ «/l“ 


sin 6 /r/A\' , r , /<fa\* 


Cela posé on aura d'abord les deux équations 

F> — FV = >(N — FM) , 
G>— G>'= >(— GM); 


dans lesquelles, si on fait /u' = k'ft, to—hfi, on aura 

,j _ rr,M 4- G*(N — FM) , _ aF’G' — aF'G" 

* f GM + G'(5f — FM) » " — F GSI + G' (N — FM)’ 

ou en faisant les substitutions et mettant les valeurs M— R cos A , 
N —R sin A, 


k'-. 


-*sinAco»AcosA. -f-sinicosA (jfj) “f- 2 tangA 

sinAcosAcosA . ^ -(-sin A. 


43. Ces valeurs supposent nr=o : on les rendra indépendantes 
de celte supposition en mettant A — a ou Ai à la place de A. 

Soit donc A — a=Ai , et pour abréger 


^ daddb — dbdda . , . dadb * . * • . s/a 3 

P = f- a taqg b . — gj- -f- sin A cos b ^ , 

Q = sin b cos b cos Ai . ^ + sin Ai . 


on aura 


co *= 4 -È-*' = 


~ dda _ . . 

— Q-^t — P si» Ai 

</n 
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Ces coefliciens donnent les valeurs de^u' et ai exprimées par le 
moyen de fi, savoir u' = A'p, a>=zhfi- il faudra ensuite déter- 
minerai et r par les équations suivantes, conformes à celles de 
l’art. 4 a : 

i* = R*-f^.aRcos 

hft — ÏÏ5, 

m' = M' + Hfi, j M' = — A i — e sin 'F cos À r , 

n> — N' + Fy, | N' = - ‘ cos Ai — -«sin 'F sin Ai , 

P' = (GA' + G>, 

w'*+ „'* + p'* = - _ i. 

* r r a 

Substituant les valeurs de ni, ri, p', dans la dernière équation, 
on aura le résultat suivant 

» 1 ‘ • * , * 

r = S + g “ 1 + + C P* > 

où l’on a fait 

B = M'A' + N'F', 

.(a) C = A'* + F'* + (GA'4-G')*; 

i 

le problème est donc réduit à la résolution des trois équation* 

■+■ hfi, 

(3) 7- = R*4- m . 2 RcosAi^-^, 

7 = l + K ~* 1 + aB ^ + C ^*' 

44- Si 1 on ne fait aucune hypothèse sur la nature de l’orbite , il 
faudra déterminer r et fi par les deux premières équations, et 
substituer les valeurs trouvées dans la troisième, ce qui fera con- 


Digitized by Google 



(64) 

naître la valeur de par laquelle on jugera s! l'orbite est ellip- 
tique ou hyperbolique selon que a sera positif ou négatif. 

Si l’on suppose d’avance que l’orbite est parabolique , ce qui 

donne i = o , on aura une équation de plus que d’inconnues. Dès- 

lors il est possible de combiner ces équations de manière à obte- 
nir la solution la plus simple, et toutes les combinaisons seront 
également bonnes, si les coefliciens différentiels qui font partie 
des données , ont toute l’exactitude nécessaire. 

Une combinaison qui se présente naturellement et qui semble 
cependaut n'avoir point encore été aperçue, réduit immédiate- 
ment le problème à la résolution de l’équation du troisième degré ; 

(4) Qî + v)(r* + ^Rco»A. + -£_) = £ — i + i*. -k Q»\ 

Ainsi dans cette méthode l’orbite peut être déterminée sans 
aucun tâtonnement, ce qui est un avantage que n’ont point 
les autres méthodes ou qu’elles n'auraient que par des élimina- 
tions très laborieuses, et qui compliqueraient extrêmement les 
calculs. 

Lorsque l’équation précédente sera satisfaite et qu’on connaîtra 
par conséquent une valeur suffisamment approchée de fi , on la 
substituera dans deux des équations (3), afin d'en tirer deux 
valeurs der; si ces deux valeurs s’accordent entre elles, on en 
conclura que la troisième est satisfaite aussi , et que le problème 
est résolu. 

Si ces deux valeurs de r ne s’accordent pas , on en conclura que 
l’hypothèse parabolique n’a pas lieu et que la quantité i qui entre 

dans la troisième des équations (3) , ne peut être négligée. Dans 
ce cas il faudra , comme nous l’avons déjà dit , résoudre les deux 
premières équations par les fausses positions , ce qui sera toujours 
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très facile, va le pea de complication de ces équations ; on con- 
naîtra ensuite par la troisième la valeur de ^ qui indiquera une 

orbite elliptique ou hyperbolique , selon qu'elle sera positive ou 
négative. 

Le cas de l’orbite elliptique peut se rencontrer dans quelques 
comètes et surtout dans de nouvelles planètes, qui, n'ayant pas 
encore été observées, seraient prises d'abord pour des comètes. 
La méthode s’appliquerait d’autant mieux à ce dernier cas, que le 
mouvement des planètes supérieures est en général beaucoup 
plus régulier que celui des comètes , et qu’ainsi les coefficiens 
différentiels pourraient être déduits des observations avec une 
précision suffisante. 

45. Examinons maintenant les différens cas qui donnent lieu à 
des remarques particulières. 

1 °. Si la quantité h est négative, on ne pourra satisfaire aux 
deux équations 

F = 5* -'*(-*)’ 

r* ==R‘ + ft( 2 R cos Ai) 

que lorsqu’on aura ( — h) < 5 En effet, si l’on considère les 

courbes décrites d'après ces équations, avec des coordonnées et r 
dont les directions soient les mêmes dans les deux courbes, il ne 
pourra y avoir d’intersection entre elles, qu'autant que la tan- 
gente au point commun, où l’on a /*=o et r=R, sera plus 

petite dans la première que dans la seconde, ce qui donne 

. i,R‘ . k - t i\ 3 cojAi 
(— A)-g-< cos Ai, ou (—/»)< — R - ~ . 

Si cette condition n’est pas remplie, la première des équa- 
tions (3) devra être rejetée , à cause de l'inexactitude du coeffi- 
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eient h, et on ne pourra avoir qu'une solution médiocrement 
approchée par les deux autres où il faudra supposer - =o. 

a°. Si la quantité h est positive, l'équation ( 4 ) aura toujours une 
racine positive, parce que le premier et le dernier terme de cette 
équation seront de signes contraires. Dans ce cas on aura néces- 
sairement r< R, et par conséquent il faudra que aRcosAi soit 

négatif et plas grand que , de sorte qu’on aura 

< a R cos* b cos ( 1 8o* + A i ). 

3 *. Si h — Oy on aura r=R et /t=aRcos*icos(i8o‘ , -f-At ). 
Dans ce même cas il faudra qu’on ait aussi i = aBju-|-C/t**. 

4°. Si h , sans cire absolument nul, est très petit , on aura à très 
peu près r=R et n— 2Rcos*£cos(i8o*-t~Ai), ce qui facilitera 
beaucoup la résolution des équations du problème. 

5 *. Si l'on a Q=o, ou sin Jcosicos Ai . ^-f-sinAi . ^ =o , il 

faudra qu’on ait en même temps P=o, sans quoi fi serait infini. 
Ce cas est celui où le lieu du soleil serait situé dans le même 
grand cercle que l’orbite apparente. Le lieu du soleil devient donc 
le nœud de l'orbite apparente à l’instant où la longitude de la 
comète est a, et sa latitude b : appelant < l’inclinaison de cette 
orbite , on aura 

lang à= lang i sin ( a — A ). 


Différenciant celte équation en faisant varier seulement a et b, oa 
en tire 

^ = sin b cos b cot (a — A).~, on Q = o; 

différenciant ultérieurement cette équation, pour achever de 
déterminer le plan de l’orbite apparente, on aura P=o, comme 

cela est nécessaire pour que l'équation RA=^ puisse subsister. 

Hais elle devient indetenkinée et ne donne point la valeur de h. 
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46. Dans ce cas par conséquent les deux équations entre /*', /* 
et a*, se réduisent à une seule, savoir FV — jN», ou 


— ï- Rsin Ai — 
r*'= Ja • 

s 

Désignons cette valeur par /»'=K.p-f-La» ; alors le problème ne 
pourra plus être résolu que dans l’hypothèse parabolique, et les 
équations réduites au nombre strictement nécessaire, seront 


r* = R* +|«(aR cos Ai) + 

i ^5 

»— i_i- K— ^ 
r» R” K ‘~ 


fji — Km + L», 
m' = M'+m'. 
n' = N' + F>, 
p' = G m' + G'fA , 
m'* n'* -f- /• = *. 


da 

3F 


» 

A i'i 


_ — JR «in Al 

33 » 

ar 


l = - 

j'i 


> V 


‘n i , 


On voit qu’en donnant une valeur à /s, on connaîtra successive- 
ment r, a, m', ni , n',p', et par ces trois dernières on aura une 
seconde valeur de r qui devra s'accorder avec la première. 

La solution est toujours possible, mais elle est plus difficile à 
obtenir daus ce cas particulier que dans le cas général ; une sem- 
blable difficulté se fera sentir, ou même deviendra plus graudo 
si les quantités P et Q, sans être loul-à-kil nulies, sont très 
petites. 

Au reste dans tous les cas, dès qu’on aura satisfait aux équa- 
tions du problème, il sera très facile d'avoir les élémcns de l’or- 
bite au moyen des quantités m, n,p, m ’ , n! , p ' , ainsi que nous 
l’avons fait voir daus tous les exemples de notre première mé- 
thode. 


I 
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47- Telles sont le* formules que l’on devra suivre dans la pra- 
tique de cette troisième méthode; elle ne laissera rien à désirer, 
toutes les fois qu’on pourra déduire avec sûreté des données de 

l’observation, les valeurs des coefficiens différentiels 

Voici au reste comment on pourra déterminer ces coef- 

ficiens d’après un certain nombre d’observations qui ne devra 
jamais être moindre que trois. 

Soit a la longitude de la comète pour une époque déterminée /, 
qui sera celle de l’observation à peu près moyenne entre toutes 
les autres; on pourra supposer pour le temps y'-f-f peu différent 
de j , la longitude représentée par la formule 

o -f- a! t H — aY - 4 — --s a!"0 -f- etc. 

' a 1 a. 3 

Appliquant cette formule à la longitude particulière a-\-ta qui 
répond au temps j + 8 , on aura l’équation de condition 

Ç = al + i 0*8 -f- ^ «'"fl* + etc. 

Dans celte équation le temps 8 , exprimé en jours moyens , devra 
être multiplié par le nombre dont le logarithme est 8.a3538ai , 
et il faudra réduire en même temps l’arc <f« en parties du rayon 
pris pour unité. 

Autant on voudra combiner d’observations, autant on aura 
d’équations semblables à la précédente, par lesquelles on déter- 
minera un pareil nombre de coefliciens «', a", «'", etc. 

La même opération servira à déterminer les coefficiens dif- 
férentiels de la latitude (*). 


(*) Il serait facile d'éviter la résolution de ces équations , car l'expression 
générale de la longitude pour un temps quelconque I , se déduit aisément 
des longitudes a, ai, oa, a3, <j<f, etc., qui répondent à des temps déter- 
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48 . Par exemple si , parmi les six observations dont noos avons 
donné le tableau, pag. 48 du i* r suppl. , on choisit les quatre des 
16, 18, a 5 et 3 o novembre, celle du a 3 étant prise pour époque, 
on aura les résultats suivans 

a = a 4 * 4 ,, 4 *» ' A = 37* a 5 ' 35 ', 

£=- [9-948a69] , "=-[ 9 . So55 7 o] , 

= -[0.675334], ^ = -[0.783466]. 

A ces données il faut joindre celles qui sont relatives au lieu 
du soleil , savoir : 

R— [9.994*071], 4 = i8o°-f-A — 99°35 ’=i4»*33'3o", 
A= 61' 8' 3 o* , M'= — [9.7856947], 

Ai=A — a= 36 . 37.36, N'= [9.9078531]; 

il eu résultera les valeurs logarithmiques suivantes 

Q = — [9.8371570], B= [o.i 5 a 7 4 oa], 

P = [o.a 3 t 6645 ], C= [1.3006704], 

h = — [o. 4 >o 5 oo 4 ], .. ; 

À' = — [o. 5446 oa 3 ]. 

Cela posé , nous aurons à résoudre les trois équations 

* » 

~= [0.0173787] — [o. 4 io 3 oo 4 ], 

r* = [9 . 988414*] + H- [° • *oo 6555 ] + /tt* [o . 1 o 356 a 1] , 

; = [ 9 - 7 I0 47 85 ] + P[°- ,52 74 o 3 ] 4 -M*[°- 99964 «> 4 ]; 


mini-s O, 8', fl', etc.; cette longitude pour le temps t est en général de 
la forme a + (qt — <j)P + (aa — a)Q •+• (a3 — o)R + etc. , P, Q, R étant 
des fonctions de t, f, l', 6”, etc., toujours connues o priori. Mais ces expres- 
sions donnent lieu à des calculs plus compliqués que ceux que nous indiquons.- 



( 7 ° y 

On ne peut , dons ce cas , faire aucun usage 4 e la première 
équation, ni par conséquent de l’équation ( 4 ), parce que ( — A), 
dont la valeur est [o. 4 io 3 oo 4 ], excède la limite 

= [0.40571 13], Ainsi ce n'est que la combinaison des 

deux dernières équations qui peut conduire à une solution ap- 
prochée. On en lire 

P- — [ 9 -i 44 i0 70 ], r= [0.0430940] ; 
calculant ensuite les quantités m, n,p , etc., d’après les formules 

ra=M + p, m' = M' tfp. , 

n = N +/u = N , n‘ = N'+F>, 

P =#* =fttang b , p' = (GA' -f- G')f* , 

on trouve les valeurs logarithmiques suivantes : 

m = [9 . 969B793] , (— m') = [o.o4°935o], (mn — m ' n ) =-+- [o. io83884], 

« = [9.768156»], n' = (J. 8357395 ], (mp— m» = — [9.41909.30], 

p = [8.85gaaoi], (—p') =[9.5640169], (np — n ' p ) = — [9.4332755], 

d’où résultent les élémens de l’orbite comme il suit : 

logn = 9.950875g. 

Passage au périh .. . déc. 29.72972, 

Inclinaison... i6'n' 8''6 * 

Lieu du Q . . . 249. 53 . 56,4 
Lieu du périh. . . 108.11.14,7. 

On voit que ces élémens diffèrent considérablement des élémens 
connus , et qu’ils ne peuvent être considérés que comme une 
première approximation assex grossière; c’est ce qu’on pouvait 
prévoir par l’inexactitude du coefficient h, qui a rendu fautive 
la première des équations ( 3 ). 

49. Pour qu’on puisse mieux juger du peu de précision avec 
laquelle les coefficiens différentiels de la longitude et de la la- 
titude peuvent être déterminés par la méthode ordinaire d inter- 
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polation , nous allons chercher , d'après les élémens corriges de 

la page 48, 1 er Suppl. , le lieu de la comète pour l'époque 
de novembre, a5. 3aa4 1 -|- 6 , 8 étant un temps très petit dont 
nous négligeons seulement les puissances supérieures à la se- 
conde, afin que du résultat contenant l’expression générale de 
la longitude et de la latitude, on puisse déduire les coeflicicus 
différentiels cherchés du premier et du second ordre. 

Époque donnée... novembre a3.3aa4' -f- 8 
Passage au périhélie... décembre 3i.a855t 

différence I = 37.96310 — fl; 

convertissant ce nombre de jours en arc du moyen mouvement 
de la Terre, et faisant j = 57 , ,<j 65 io= o. 6530460, on anra 

Voici maintenant les diverses parties du calcul que nous de- 
vons exécuter. 

I. Soit 4 l'anomalie qui répond au temps t , et n=tang;4 j 
on aura Icquation 

v , (an 3 )-(« + i^, 

d'où résulte (!+«*) ^=—15, ou g = ; ensuite, 

— — atangi-4-.y. Soit ensuite i=/; ou aura par b table 
des comètes, 4 = 53* 36' 54", et en même temps, 

§ = — [o.oa 7 5 9 5a], jr = — [0.0597703]. 

Donc, si on appelle 4* l’anomalie qui répond au temps Isa J 0, 

on aura 

4 ' = + “ 6 [o.o 0 a 7 5()5 a] — 8* [9.7587403]. 

II. L’équation donne Iogr= 0.0490374; ensuite 
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on aura par la différentiation, 

“s:=.i.. s ;+.?+ 

êt en substituant les valeurs trouvées = [9.3689294] ; 

= [ 8 . 76 . 1853 ]. 

Donc, si l’on appelle r' le rayon vecteur qui répond au 
temps l — j — 0,1 aura 

log t 1 = 0.0490374 — 6 [9. 3689294] + fl* [8. 4601 553 ]. 

111 . Connaissant la distance du périhélie au nœud descendant 
T = 38 " 5 o' 5 ", et l’anomalie 4 = 53 * 36 ' 54 ", on aura l'argu- 
ment de latitude 

4— r = . 4- 46 ' 49 "- 

D’après cette valeur de 4 — r et celle de l’inclinaison 
1= 16* 3i' 10", on déterminerais latitude hcliocentrique € et 
la longitude a, par les équations sin Ç = sin I sin (4 — T), 
tang (Ni — a) = cos I lang (4 — r] ; on trouvera ainsi 


sinS = [8.8605716] , 
cos 6 = [9.9988544] , 


Ni — a = i4* n'43"3 
Ni = 70.33.54 

a = 56.ai .5o,7 
A = 61. 8.3o 


A— a = S = 4.46.39,3. 

De l’équation sin Ç= sin I sin (4 — r), on déduit 
= m col (4-T).ÿ =[o. 3440536] , 

= "» «t (4— r) . -^=-[0.9761717]. 
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Donc , si fi' est la valeur de fi qui re'pond au temps te=/ — 0 , 
ou aura 

Zsinfi' = 8.8605716 — 0 [o.a 44 o 336 ] — fi* [0.6751417] » 
d'où résulte 

Zcosfi' = 9.9988544 +0[7 -9674680] — 6* [ 8 . ioi5573]. 


Ajoutant la valeur connue de log r 1 , on aura 

L (/ sin fi') = 8. 9095990 — 8 [0.398386a] — fl* [0.6734864] , 
Z(r'c os fi') = o . 04788 1 8 — 0 [7 . 35 1 3467] -f- 8* [8 . 3097478]. 

IV. De l’équation tang (Ni — *) = cos I laDg (4 — T) on 
déduit 

g^sin(a+— ar)=r— ^ «in(aN — a«) , 

~sia(a 4 — ar)r= — «io(aN— a«) + ^^'[co»(aN— a«)— co<!h| — ar)J 
ce qui donne 

~ = — [o.on 58 oo], ^ = [0.0261734]. 


Donc, en appelant *' la longitude pour le temps t—j — 6 , 
on aura 

#' =: 56 * ai' 5 o ",7 + 8 [o.on 58 oo] + 0 * [9.7351434]; 


on a semblablement pour le lieu correspondant de la Terre, vu 
du Soleil, 


A' = A — 8.~ ; 


. «MA 
■~i F» 


or, par les formules de la page 10 du premier Mémoire, on 
trouve 


dk t/(*“ »*) 1 ddk VO — •*) dR i»inf- 

dt — R* • dr ~~~~ RÎ • dt ~ 


R 1 ». 
10 
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substituant les valeurs t = [8.3a5o5i] , T = i 4 »° 53 ' 3 o", ou 
aura 

A' = 6i° 8' 3 o" + 6 [o.ou 5 a 46 ] + 0 * [ 8 .o 36 oaa 8 ], 

. a.' = 56 . ai . 5 o ,7 + 0 [o. oi i 58 oo] -f- 0 ’ [9.7351424]» 

A' — *'=S' = 4*46-39,3 — 0 [6. 1173500] — 6* [9.7161G56] j 

il en resuite 

logsinS' = 8 .gao 58 a 3 — 0 [6.8339404] — 0 * [o. 43 i 856 i] , 
logcosS' = 9.9984884 H- 0(4.6771383] — 0 * [8.3760457]. 

Connaissant déjà, pour le temps y — 0 , la valeur de SQ'=r'cos£', 
on en déduira celle SI' = SQ'cosS' et celle de QT=SQ'siuS', 
dont les logarithmes sont 

logQT = 8.9684641 — 0 [9.3536611] 4- 0 * [0.4992428], 
log SI' = 0.0465702 — 0 [g. 35 i557d] -f- 6* [8 . 5463820]. 

De ce dernier on déduit, en passant aux nombres, 

SI' = 1.1 126797 — fl [ 9 - 7 5 99 a 34 ] + G*[ 9 - 5 7779 58 ]. 

On a encore, par les formules citées, R' ou 

ST' = 0.9867500 — 0^8.0187053] -f- 0 ’ [7.8295025]; 

, donc, 

TT = 0.1359297 — ô r g. 7519705] + B*[ 9 . 36553 o 5 ], 
log TT = 9.1001282 — 0 [0.2896264] — 6’ [0.5536589]. 

Enfin en appelant a la longitude qui répond au temps y — 0 , 
on a tang(A' — a!) = 2 tj. ; donc 

L tang (A' — a’) = g . 868535 <)+â[o ■ a 56 a 65 i ]-M*[9 . 94^91 81], 

A' — a' = 36 “ 36' 4a , '+'i[o.2777gi5]-|-0 , [o.3i76i62] , 

A' = 61. 8 . 5 o -f- 0 [o.oi i5a46]+0*[8.o36o228] , 

«' = 24.41.48 9389714]— 6 *[o. 3 i 535 g 4 ]. 
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Celte dernière valeur donne les deux coefficiens différentiels 
de la longitude, savoir, 

^==—[9.9389714], ^'=—[0.6165694]. 

Pour déterminer la latitude , il faut faire le calcul de la for- 
mule tang b' = cos (A' — a'); on en déduira 

L tang b' = 9.7149577— 9 [9.8115034] — fi’ [0.4781 180] , 
b 1 = 27” a 5 ' 4"8 — 8 [9.7849805] — 8* [0.4096529] , 

et enfin les coefficiens différentiels de la latitude, 

^= — [9:7849805], ^= — [0.7106829]. 

On avait trouvé, par l'interpolation de quatre observations les 
plus proches de l'instant pris pour époque, 

£ = - [9-948369], ^ - [o. 673534] , 

= - [9-8o55 7 o}, ^ - [0.782466]. 

Ainsi on voit que l’interpolation ne donne pas avec assez 
de précision les coefficiens différentiels cherchés, surtout ceux 
du second ordre, ce qui nuit beaucoup au succès de la mé- 
thode. 

\ 

5 o. 11 n’y a pas de doute qu'avec les coefficiens différentiels 
que nons venons de trouver, l'application de la méthode ne 
fasse connaître la véritable orbite d'une manière suffisamment 
approchée. C’est ce que nous allons vérifier par le calcul suivant. 

11 faut d’abord remarquer que les élémens corrigés sur les- 
quels nous avons établi les calculs précédées, donnent le lieu 
du a 3 novembre un peu différent du lieu observé; en effet, 
nous avons trouvé par ces élémens la longitude a — 24° 4 ,, 48 ’ , , 
et la latitude b = 27° 2 5 ' 4", 8, tandis que le lieu observé donne 
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o=3 4 * 4 ,( 4 ", et b = 37* 35' 55 ", ce qui donne à peu près les 
memes erreurs qu’on a notées dans le tableau de la page 47- 
Pour l’application de la méthode, il faudra donc supposer 
a = 34* 4»' 48", * = 37*35' 4", 8> on a d’ailleurs pour le lieu 
du Soleil les données 

R = [9.0943071], A = 6i* 8' 3 o", 

* = 141* 53 ' 3 o", Ai=A— o = 36 .a 6 . 4 a; 

substituant ces valeurs, ainsi que celles de s = [ 8 .aa 5 o 5 i} , 
dans les formules 

M = RcosAi, M' = — ^ '-~*' ^sinAt 
N = RsiuAi, N' = Ç -• ■ ‘ -)cosAi 
on en tire 

M = [9.8996940] , M' = - [9.78557.8 6], 

N = [9.7680307], N' = [9.9079030 5 ]. 

On trouve ensuite les valeurs logarithmiques des quantités qui 
servent à former les équations du problème, comme il suit : 

F ' = S=-C9-95897«4L A=-[o.53.6356], 

G = tangè= [ 9-7 >49677] » *'=—[0.48943.8], 

G ~d^’ 2 i=-fe- 888 47 6 7 J> G*+G'= — [0 . 375576a 5 ], 
P= [0.1373743], B= [0.073.754], 

Q= — [9.81.43.5], C= [..3018967]. 

Cela posé, nous aurons à résoudre les trois équations 

4 j= [0.0.73787] — fi [o.Sa. 6356 ], 

r* = [9. 9884.43] + M [o- 200734°] + ^*[0.1034963], 

f = b-7 10 479 5 ] -1- H- [0.072.754] ■+■ fi * [0.9008667]. 


— «sin'PcosAr, 

— isin'Psin Ai , 
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Si. Ces trois équations s'accordent à donner à très peu près 
les mêmes valeurs de r et de p. Par les deux dernières, qui 
doivent être les plus exactes, on trouve 

f* = [g-’Ç^Sgaa], r = [0.0490209]. 

Pour déduire de cette solution les élémens de l’orbite, on 
calculera d’abord les valeurs logarithmiques des quantités m , 
n, p , m ' , n', p' , ce qui donnera les résultats suivans : 

m=[g. 977859a 5 ], m’= — [0.038789a a], mn ' — m'n= [0,1073459 a 3 , 
n= [9. 7880307], n = [9 837904a], mp — m'p= — [9.4330955 7J, 

P=C8. 9095499], /=— [9.5701684 5 ], np — n ’ p — — Q). 43538 a 5 3 ]. 

Cela posé, la distance périhélie IT se trouve directement par la 
formule 

an = {mn!— m'/i)* + ( mp m'p )• -4- {np’ — n'p) % , 

d’où résulte n = [9.9502708 4] ; l’inclinaison I, déduite de 
la formule cos 1 = donne cos 1 = [9.9816965] , et 

I = 16° 3 i' 4 ", 8 . Le lieu du nœud S est donné par la formule 
tang(S — a) = = [0.0130869 ®]» d’où résulte 

S — a = 45 ’ 5 i' 4 7 "5 

a = 34.41.48 

S == 70* 33 ' 55''3 
ou S == a 5 o. 33 . 55 , 3 . 

La seconde valeur est celle qu’il convient d’adopter, parce 
que cos(S — a) doit être de même signe que mp' — m'p. 

De l’équation n=rcos*î 4 , on déduira C0S;4 = [9-95oGa497], 
W= a6 ° 48' 29", 5 a , 4 = 53 “ 36 ' 4 7 '',o 4 . 

A cette anomalie répond, dans la Table des comètes, le 
temps T = 45 .o 7 4 g 58 =[i. 6539353 ], d’où l'on déduit le 



1 


C 7® ) 

temps < employé à décrire l’anomalie 4 > 

t = n*T = 37.961344 

Epoque du lieu observé. . . nov. a 5 . 82241 
Passage au périhélie... déc. 51.283754. 

Enfin l’cquation tang 9 [9.7917145] donne la longi- 

tude héiiocentriqnc <p = 3 i* 4 °' 4' , »4 » d’où l'on déduit aisément 
Je lieu du périhélie = 109* a 3 ' 34 ",!. 

En rassemblant tous ces résultats , on aura les élémens de 
l’orbite comme il suit : 


9.9502708 
déc. 31.28375 

16° 3i' 4"8 
a5o.55.35,5 

109.23.34,1 

On voit que ces élémens s’accordent presque entièrement 
avec les élémens corrigés de la page 48, i cr suppl., les diffé- 
rences étant très légères et devant être attribuées, soit à la lon- 
gueur des calculs, soit à l'imperfection des données qui ont 
servi à déterminer les coefliciens M' et N'. 

Paris, le 10 août i 8 ao. 


Log n 

Passage au périhélie 

Inclinaison 

Lieu du nœud asc.. 
Lieu du périhélie. . . 
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NOTE par M. ***. 


Au commencement de la dissertation précédente, page 4 , l’auteur parlant 
de sa Méthode des moindres cariés, s'est borné à rappeler qu'il l'avait pu- 
bliée, pour la première fois, en » 8 o 5 , à la (in d’un mémoire sur les comète* , 
qui porte la da'e du mois de mars de cette année. Cependant, comme un 
géomètre fort célèbre n’a pas hésité à s’approprier cette méthode , dans un 
ouvrage imprimé en 1809, nous croyons devoir nous arrêter un instant sur 
cette prétention que tout lecteur impartial pourra qualifier ensuite du nom 
qu’il jugera convenable. Voici comment elle est énoncée : Après avoir exposé 
cette méthode ou ce principe comme lui appartenant, ce géomètre ajoute: 
u Ceterum prinùpium nostruni quo juin indè ab anno 1795 usi sumus, nuper 
« etiam à clar. Legendre in opéré Nouvelles Méthodes, etc prolatum est , etc. « 
Si cela n’est pas très décisif, c’est du moins fort clair, et surtout f«rt com- 
mode. Dans ce système, l’histoire des sciences s’écrira beaucoup plus aisé- 
ment; une découverte pourra bien ne plu? appartenir à celui qui l'aura faite; 
mais qu’importe I Elle appartiendra toujours à quelqu'un, à celui qui aura 
trouvé bon do la revendiquer sam titre , à partir de l’époque la plus reculée. 
De bonne foi, un pareil système est-il admissible ? On ne l'avait pensé jus- 
qu'ici dans aucune circonstance , comme l’Histoire des Mathématiques en fait 
foi à chaque page; on regardait la propriété d’une découverte comme inva- 
riablement assurée à celui qui, le premier, la mettait au jour, et toute prétention 
contraire à ses droits, pouvant l'txposer lui- même à des soupçons d’une nature 
injurieuse , exigeait l’appui de titres précis et authentiques. Au-si dirions-nous vo- 
lontiers que celle dont nous parlons était, en i8cg, d'une espèce tout-à-fait 
nouvelle , si nous ne trouvions qu’en 1801 le même géomètre s'essayait dès-lors, 
dans ce genre, d’une manière il est vrai bien plus imparfaite. On avait vu paraître 
dans les Mém. de l Ac. des Sc. pour 1785, un théorème ( la loi de réciprocité ) 
dont les conséquences fesaient faire des progrès réels à l’analyse des nombres, 
sans que l'inventeur, on en convient, l'eût démontré rigoureusement dans tous 
les cas. Or , comment le géomètre dont nous parlons s’exprimait-il, à ce sujet, 
teize ans après? 11 établissait d’abord ce théorème, par lui nom me fonda- 
mental à cause de son importance, puis il ajoutait: a Theoiema fundamen- 
n taie quod sane inter e/eganiissima in hoc genere est referendum , in câdem 
t> forma in quâ suprà proposilum est , a nemine hue usque fuit prolatum, rt 
Ou trouve cependant deux pages plus loin : a Clar. Legendre in egregid tracta - 
tione Ac. des Sc. 1785, pervenit ad theorema quod si rem ipsarn spectas cum 
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i» theoremate fundamentali idem est. p Voilà un aveu qui ne s'arrange guère 
avec ce qui précède, et un bien singulier théorème qui, publié en 1785, n'avait 
pourtant été donné par personne en 1801 ! Mais du moins à cette époque on 
ne disait pas encore , en se donnant dix ans d'avance , comme on a fait 
depuis en 1809 ; u theorema nostrum quo jam inde ab anno 1775 nsi sumus, n 
Quoi qu'il en soit , au surplus , de ce système sur la propriété littéraire , il 
ne sera peut-être pas inutile de rappeler en finissant, que la fameuse série de 
Lagrange , pour le retour des suites , quoique non démontrée par son au- 
teur , n’a pas cessé de porter le nom de Théorème de Lagrange , et que nul 
jusqu'ici n'a songé à lui en disputer U propriété. 


De l'Imprimerie de HUZARD-COÜRCIER , rue du Jardinet , n® ta. 
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'ERRATA du second Supplément. 


Page 


3, 

ligne 

8, 

-, lisez - 

a a 

s 

, m 

4 , 


6, 

Sy, lisez, S“T 


7» 


6, 

*/<...£, lisez Ac. . 

/ 
" C 

>3, 

ligne 

antépénultième, <aA\ t 

lisez cos A i 

>5, 

ligne 

20, 

deux ans , /w« deux cas 

>7. 


5, 

négative , lisez de signe contraire à h 

ai, 


10, 

■jR cos ( 1 8o° -f- A ■ 

— n) , lisez a R coj* b cos ( 1 8o* -f- A -—ai) 

a3. 


18, 

lisez log a' = log ^ 

*' \ . « , 
5) -f- 6V Û+T * 

a5. 


«4, 

t ang (l±i _ jy\ 

, //ses tang T” - ~~ “* ^ 

07, 

ligne 

pénultième , â et 6', lisez A et a' 

36, 

ligne 

9, 

cosi4— 9 - 98 , 

, lisez cos £ 4* = g . 89 . . . 

63, 



La troisième des équation} (3) doit être lue ainsi : 




- = i -f- — — 
r a W{ 

• 1 -f- tsBfâ -f- Cft* 



3, 

à la place de la première phrase, lisez , 


fl®. si la quantité h est positive, il faudra, conformément à l'art. ia, 
que la quantité cos Ai + 3 A/î» soit négative. Dans ce cas, etc. 

J 2 , lisez a + fl , l+ i aV -t- ^73 a"* 3 -f“ etc. 
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